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ABSTRACT 


The  aqueous  solution  chemistry  of  trimethyllead ( IV) 
and  of  the  trimethyllead (IV)  complexes  of  selected 
organic  ligands  has  been  investigated  by  proton  magnetic 
resonance  spectroscopy.  The  formation  constants  for 
the  (CH-^^Pb"1"  complexes  of  hydroxide  ion  and  of  the 
various  sulfhydryl  and  carboxyl  functional  groups  of 
2-mercaptoethanol ,  glycine,  N-acetyl-D ,L-penicillamine , 

D , L-penicillamine ,  L-cysteine  and  glutathione  were  deter- 
mined  from  the  pH  dependence  of  the  chemical  shift  of 
the  methyl  protons  of  trimethyllead.  One-to-one 
stoichiometry  was  dominant  in  all  experimental  cases. 

In  addition,  weak  two-to-one  metal-to-ligand  complexes 
were  observed  for  the  trimethyllead  complexes  of  hydroxide 
ion  and  of  the  sulfhydryl  groups  of  2-mercaptoethanol 
and  glutathione.  No  evidence  of  amino  group  complexation 
by  trimethyllead  was  found.  The  formation  constants  of 
the  carboxylate  complexes  and  of  the  various  two-to-one 
metal-to-ligand  complexes  are  small,  while  those  of  the 
sulfhydryl  complexes  of  one-to-one  stoichiometry  are 
generally  several  orders  of  magnitude  higher.  The  extent 
of  complex  formation  is  strongly  pH  dependent  due  to  the 
competitive  effects  of  ligand  functional  group  protonation 
at  neutral  or  low  pH  and  (CH^^PbOH  formation  at  high  pH. 
This  pH  dependence  is  illustrated  by  calculations  of  the 
species  distributions  and  the  conditional  formation 
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constants  of  the  various  complexes  as  functions  of  pH. 
The  trimethyllead  chemical  shifts  of  the  sulfhydryl 
and  carboxylate  complexes  are  found  further  upfield  and 
the  formation  constants  increase  as  the  pKa  of  the 
ligand  functional  group  increases.  To  some  extent, 
other  factors  such  as  steric  and  electrostatic  effects 
also  appear  to  influence  these  parameters. 
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CHAPTER  I 


INTRODUCTION 

Inorganic  lead,  in  the  forms  of  metallic  lead  and 
lead  salts,  has  been  familiar  to  mankind  since  ancient 
times.  In  contrast,  compounds  which  include  a  lead- 
carbon  bond  were  unknown  until  1853,  when  Loewig  pub¬ 
lished  a  series  of  papers  describing  the  synthesis  of 
several  organolead  compounds  (1,2,3).  Included  were  the 
first  reports  of  organolead  salts  of  the  type  R^PbX, 

X  being  an  inorganic  anion.  In  1860,  Klippel  prepared 
a  number  of  similar  compounds  using  organic  anions  (4) . 
From  these  origins,  organolead  chemistry  has  developed 
into  one  of  the  largest  areas  of  organometallic  chemistry, 
exceeded  only  by  those  of  tin,  silicon  and  mercury  (6) . 
This  prominence  is  primarily  due  to  the  utility  (and 
hence,  large  scale  production)  of  various  organolead 
compounds  as  antiknock  additives  for  internal  combustion 
engine  fuels  (7) ,  a  property  first  reported  in  1922  (8)  . 

The  most  recent  comprehensive  reviews  of  organolead 
chemistry  (5,9)  and  the  annual  reviews  published  by  the 
Journal  of  Organometallic  Chemistry  demonstrate  the 
preponderance  of  synthetic  and  non-aqueous  chemistry  in 
the  organolead  literature.  Relatively  little  is  under¬ 
stood  of  the  aqueous  coordination  chemistry  of  organolead 
compounds;  this  is  especially  true  of  the  cationic 


1 


triorganolead  species.  However,  recent  studies  of  the 
environmental,  toxicological  and  chemical  behaviour  of 
such  compounds  have  emphasized  the  need  for  a  greater 
understanding  of  the  properties  of  aqueous  triorganolead. 
Partial  fulfillment  of  this  need  is  presented  in  this 
thesis,  a  study  of  the  aqueous  solution  chemistry  of 
trimethyllead  (the  simplest  of  the  triorganolead  cations) 
and  of  its  complexes  with  a  number  of  organic  ligands. 

A.  Sources  and  Toxic  Action  of  Trimethyllead 
1 .  Trimethyllead  in  the  Environment 
Several  independent  processes  are  thought  to  be 
responsible  for  the  environmental  formation  and  accumu¬ 
lation  of  triorganolead  compounds.  Atmospheric  contamina¬ 
tion  with  organolead  compounds  arises  primarily  from  the 
vaporization  of  unburned  gasoline  additives  (10,11), 
the  main  sources  being  open-air  transfer  of  fuel  between 
containers  and  incomplete  combustion  of  fuel  and  waste  oil. 
In  addition,  organolead  compounds  are  being  employed  as 
biocidal  agents,  as  chemical  catalysts  and  as  stabilizers 
for  PVC  plastics,  as  well  as  in  lubricating  oils,  in 
electroplating  processes  and  in  a  number  of  polymerization 
procedures  (10,12).  The  manufacture,  use  and  ultimate 
disposal  of  such  products  results  in  further  contamination 
of  the  atmosphere,  soils  and  natural  waters. 

Triorganolead  compounds  form  only  a  portion  of  the 


total  organolead  released  to  the  environment.  However, 
some  of  the  other  organolead  compounds  are  thermodynamic¬ 
ally  unstable  (13) ,  being  decomposed  by  light  or  by  air 
oxidation  at  fairly  low  temperatures  (14) .  Toxic  tri- 
organolead  compounds  are  rather  stable  intermediates  in 
such  decompositions. 

It  should  also  be  noted  that  organolead  compounds 
account  for  only  a  fraction  of  all  lead  pollution. 
Extensive  commercial  use  of  lead  in  many  forms  (10)  has 
resulted  in  significant  environmental  concentrations  of 
inorganic  lead.  A  major  contributing  factor  is  the 
complete  combustion  of  organolead  gasoline  additives, 
with  the  resultant  formation  of  inorganic  lead  compounds 
(12) .  Indeed,  only  a  few  percent  of  the  total  lead  in 
urban  atmospheres  is  organolead;  the  remainder  is 
inorganic  lead  particulate  matter  (15) .  However,  many 
of  these  lead  compounds  can  be  methylated  by  mixed 
microbial  communities  in  lake  sediments  (16-18) ;  a  general 
mechanism  for  lead  biomethylations  has  been  proposed 
(19) .  Since  most  of  the  lead  released  by  man  eventually 
is  trapped  in  soils  and  sediments  (20)  where  such  microbes 
flourish,  appreciable  amounts  of  trimethyllead  and 
tetrame thy Head  may  be  formed  in  this  way.  In  addition, 
tri-  and  tetraorganolead  compounds  may  be  produced  by 
disproportionation  of  diorganoleads  (18). 

The  multiplicity  of  interconversion  mechanisms  of 
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organic  and  inorganic  lead  compounds  has  caused  significant 
uncertainty  in  the  estimates  of  environmental  triorgano- 
lead  concentrations.  Nonetheless,  these  compounds  are 
sufficiently  toxic  to  warrant  concern,  especially  since 
higher  concentrations  are  expected  in  densely  populated 
areas  (10,11,20) . 

2 .  Trimethyllead  Toxicology 

Inorganic  lead  poisoning  is  relatively  common  and 
its  toxic  effects  and  treatment  with  chelating  agents 
have  been  well  documented  (11) .  In  contrast,  the  inci¬ 
dence  of  acute  organolead  poisoning  is  low;  hence,  much 
less  is  understood  of  its  toxicology  and  treatment. 

It  is  known,  however,  that  standard  chelation  therapy 
with  2 , 3-dimercaptopropanol  (BAL) ,  EDTA  or  penicillamine 
is  ineffective  (21-24)  except  for  cases  of  diorganolead 
poisoning,  which  respond  to  treatment  with  BAL  (23,24). 
Almost  nothing  is  known  of  the  effects  of  long  term 
low  level  exposure  to  organoleads  (11,21). 

Tetra-  and  triorganolead  compounds  induce  distinctive 
toxic  effects  unrelated  to  the  toxicity  of  inorganic 
lead  (22) .  Tetraorganoleads  appear  to  be  relatively 
inert  in  themselves,  but  considerable  evidence  indicates 
that  "in  vivo"  decomposition  of  these  compounds  occurs 
(primarily  in  the  liver) ,  forming  triorganolead  salts, 
the  ultimate  toxic  agents  (21-29)  .  Hence,  absorption  of 


tetraorganolead  compounds  through  ingestion,  inhalation 
or  skin  contact  will  result  in  symptoms  identical  to 
those  of  triorganolead  poisoning. 

The  toxic  effects  of  triorganoleads  appear  to  vary 
according  to  the  degree  of  lipophilicity  of  the  organo- 
lead  cation  and  of  the  neutral  complex  (25) .  It  is 
thought  that  the  primary  effect  of  the  relatively  hydro¬ 
philic  trimethyllead  species  is  to  promote  increased 
exchange  of  mitochondrial  hydroxide  ions  for  extra- 
mitochondrial  chloride  ions,  resulting  in  swelling  and 
altered  organelle  function  (25,30).  Other  documented 
effects  include  inhibition  of  mitochondrial  oxidative 
phosphorylation  (that  is,  synthesis  of  ATP),  inhibition 
of  ATP-coupled  enzymes  which  transport  sodium  and  calcium 
ions,  inhibition  of  glutathione  transferases  and  alter¬ 
ation  of  mitochondrial  membrane  permeability  (25,30,31). 
However,  the  relation  of  these  effects  to  the  dominant 
action  of  triorganoleads  on  the  central  nervous  system 
is  not  well  understood  (25) . 

All  of  the  above  effects  are  observed  primarily 
during  acute  triorganolead  intoxication,  but  other 
physiological  disturbances  may  result  from  much  smaller 
doses.  For  example,  trimethyllead  is  known  to  bind  to 
the  phosphate  backbone  of  DNA  (32) .  It  is  therefore  not 
surprising  that  trialkylleads  are  somewhat  mutagenic  in 
certain  plants  and  in  larval  fruit  flies  (33) .  No 
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teratogenicity  has  been  found  for  such  compounds  (34), 
but  studies  confirming  this  have  not  been  performed. 

In  addition,  low  doses  of  tetraethyllead  are  carcinogenic 
towards  mice  (35) ;  paradoxically,  lethally  high  doses  were 
reported  to  have  a  regressive  effect  on  mouse  tumors 
(36) .  In  man,  it  has  been  determined  that  long  term 
contact  with  tetraethyllead  inhibits  spermatogenesis  (37) 
and  may  also  result  in  profound  changes  in  the  hemopoietic 
system  (38) . 

To  date,  no  efficacious  therapeutic  agent  for  tri- 
organolead  poisoning  has  been  demonstrated  conclusively. 

(An  early  report  claiming  dramatic  improvement  with  a 
mixture  of  sodium  trithiolactate  and  sodium  thiolactate 
(39)  remains  unconfirmed.)  Fortunately,  organolead 
poisoning  is  relatively  self -limiting  in  duration: 
continuous  supportive  care  allows  eventual  full  recovery 
(21) .  The  development  of  a  rapid  and  effective  therapy 
is  nonetheless  an  important  end,  but  the  prerequisite  means 
to  this  are  a  deeper  understanding  of  both  the  toxic 
mechanisms  and  the  coordination  chemistry  of  triorganolead 
compounds . 
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B.  The  Coordination  Chemistry  of  Trimethyllead 

1 .  The  Structure  and  Bonding  of  Trimethyllead 

Complexes 

It  is  well  established  that  tetramethyllead  and  other 
R^Pb  compounds  are  tetrahedral  (40,41).  However,  cleavage 
of  one  lead-carbon  bond  produces  the  trimethyllead  cation, 
(CH^J^Pb4,  which  may  assume  a  variety  of  structures  depen¬ 
ding  on  the  number  and  kind  of  Lewis  base  molecules  co¬ 
ordinated  to  the  lead  atom.  When  only  one  base  is  bonded 
to  each  cation,  as  occurs  for  solutions  of  trimethyllead 
salts  in  non-coordinating  solvents,  the  complexes  are 
four-coordinate  pseudotetrahedral  monomers  (42-44) . 

In  most  cases,  however,  the  favored  geometry  of 
trialkylleads  is  five-coordinate  trigonal  bipyramidal  (tbp) . 
Indeed,  solid  state  trimethyllead  complexes  of  one  to 
one  stoichiometry  commonly  exist  as  polymeric  chains  of 
tbp  units,  with  the  anions  bridging  between  planar  tri¬ 
methyllead  moieties.  Examples  include  the  trimethyllead 
complexes  with  halides  (42,45),  methoxide  (45),  formate, 
acetate,  monochloroacetate  and  imidazole  (43,46),  cyanide 
(48),  thiocyanate  (44)  and  several  other  donor  species 
(49-52) .  Exceptions  to  this  rule,  such  as  the  monomeric 
thio-  and  dithioacetate  complexes  (43)  ,  are  very  uncommon 
in  the  literature.  X-ray  crystallographic  studies  have 
confirmed  the  tbp  polymeric  structure  of  solid  trimethyllead 


acetate  (53) . 


Trigonal  bipyramidal  symmetry  is  also  favored  in 
solutions  of  these  complexes  in  strong  donor  solvents. 

In  such  solvents,  the  planar  trimethyllead  moiety  is 
axially  coordinated  by  two  donor  molecules,  at  least  one 
of  which  is  a  molecule  of  solvent.  A  number  of  such 
coordination  compounds  of  trimethyllead  have  been  des¬ 
cribed  (44,49,42,54-57)  . 

The  strength  of  the  electron  donors  is  a  factor  in 
complexations  of  this  sort.  For  example,  the  work  of 
Shier  and  Drago  (57)  on  trimethyllead  perchlorate  in  pyri¬ 
dine  solution  showed  that  the  perchlorate  anion  is  not 
coordinated  directly  to  the  lead  atom  as  it  is  in  the 
solid  state,  but  is  replaced  in  the  inner  coordination 
sphere  by  a  second  pyridine  molecule.  Hence,  a  strong 
donor  solvent  molecule  may  displace  the  anion  of  a  solid 
complex  to  the  outer  coordination  sphere  during  solvation 
(58) .  Water  is  another  solvent  which  may  behave  in  this 
manner  (52)  . 

Less  polar,  less  strongly  coordinating  solvents  also 
tend  to  preserve  tbp  symmetry  in  trimethyllead  complexes, 
but  by  coordination  of  only  one  solvent  molecule  to  the 
lead  atom  (47) .  Here  the  anion  remains  in  the  inner 
sphere  of  the  complex.  Little  or  no  solvent  coordination 
occurs  in  the  non-polar  solvents,  resulting  in  a  tetra¬ 
hedral  trimethyllead-ligand  complex  (52).  However, 
perfect  tbp  and  tetrahedral  conformations  are  observed 
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only  for  certain  complexes  in  the  strongest  and  weakest 
donor  solvents,  respectively.  Intermediate  distorted 
symmetries  may  also  occur,  depending  on  the  differences 
in  donor  properties  of  the  ligand  and  the  solvent  (59,60). 

Water  is  an  almost  ideal  polar  solvent  for  trimethyl- 
lead  cations,  as  it  has  the  necessary  qualities  of  good 
donor  characteristics  and  a  high  dielectric  constant  (13) . 
In  aqueous  solution  the  (CH^)  (3IV)  moiety  is  therefore 

planar  (61)  and  is  assumed  to  be  axially  coordinated  by 
two  water  molecules  (62).  One  of  these  solvating  water 
molecules  may  be  displaced  by  certain  ligands  (58,62,63), 
often  without  altering  the  symmetry  significantly.  (For 
simplicity,  any  such  waters  of  solvation  are  omitted  from 
the  molecular  formulae  of  the  various  aqueous  trimethyllead 
species  throughout  this  thesis.)  A  consequence  of  this 
model  is  that  chelate  formation  with  trimethyllead  is 
not  favored;  this  is  consistent  with  the  ineffectiveness 
of  chelation  therapy  noted  earlier. 

The  symmetry  of  complexed  trimethyllead  is  a  function 
of  the  types  and  strengths  of  the  various  metal-ligand 
bonds.  The  methyl  groups  of  the  aquated  trimethyllead 
cation  are  coplanar  because  the  lead  s-orbital  tends  to 
concentrate  in  the  bonds  to  these,  the  least  electro¬ 
negative  substituents  (57) .  The  lead-carbon  bonds  are 

2  . 

therefore  essentially  sp  hybridized,  with  33%  s-character 
in  each  bond  (55) .  The  axial  bonds  to  the  water  molecules 
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are  formed  from  the  lead  p  and  d  2  orbitals,  which  possess 

z  z 

the  appropriate  orientation  (60) .  The  absence  of  s- 

character  from  the  axial  bonding  orbitals  results  in 

longer,  weaker  and  more  polarized  bonds  to  the  water 

molecules  than  to  the  methyl  groups  (55)  . 

There  are  thus  two  distinct  types  of  metal  bonding 

in  tbp  complexes  of  trimethyllead :  the  kinetically  stable, 

2 

sp  hybridized  and  relatively  covalent  lead-carbon  bonds 

and  the  labile,  relatively  polar  and  p  +  d  2  hybridized 

bonds  to  the  electronegative  donor  molecules  (9,13,52). 

However,  the  nature  and  hence  the  symmetry  of  these  bonds 

is  dependent  on  the  characteristics  of  the  ligands.  The 

bonds  to  "hard"  bases  (as  defined  by  Pearson  (64))  have 

a  high  degree  of  ionic  character,  whereas  those  to  "soft" 

bases  are  more  covalent  (13,52).  Distortion  of  tbp 

symmetry  towards  tetrahedral  arises  from  strongly  covalent 

ligand  bonding,  which  requires  transfer  of  s-character 

away  from  the  lead-carbon  bonds  to  the  lead-ligand  bond 

3 

in  question.  The  lead  orbitals  become  increasingly  sp 
hybridized  and  d  orbital  participation  is  concurrently 
reduced.  In  addition,  the  symmetry  may  be  affected  by 
possible  pTT— dn  and  hypercon jugative  interactions  which 
have  been  reported  for  trimethyllead  complexes  with  sulfur 
functional  groups  (65—67)  and  unsaturated  systems  (68,69) , 
although  such  effects  are  negligible  with  relatively 
"hard"  Lewis  base  complexes  (70,71) .  Steric  and  lipo- 
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philic  interactions  of  ligand  substituents  with  the  tri- 
methyllead  methyl  groups  may  also  affect  the  symmetry  as 
well  as  the  strength  of  ligand-metal  bonding.  (Such 
properties  come  into  prominence  in  the  biochemistry  of 
trimethyllead ;  selective  binding  to  macromolecules  is 
thought  to  depend  in  part  on  the  presence  of  special 
tertiary  structures  (25).) 

For  the  purposes  of  this  thesis,  the  symmetry  and  co¬ 
ordination  chemistry  of  trimethyllead  can  be  summarized 
rather  simply.  The  aquated  trimethyllead  cation  is  tbp 
in  symmetry,  with  one  water  molecule  coordinated  in  each 
of  the  axial  positions.  Complex  formation  involves  the 
displacement  of  at  least  one  of  these  water  molecules  by 
another  ligand.  Depending  on  the  donor,  the  tbp  symmetry 
may  be  distorted  to  some  extent,  but  the  central  lead 
atom  will  almost  always  remain  pentacoordinated ,  with  the 
labile  ligands  positioned  axially. 

2 .  Formation  Constants  of  Trimethyllead  Complexes 
in  Aqueous  Solution 

Quantification  of  the  binding  of  aqueous  trimethyllead 
by  various  ligands  has  been  the  subject  of  very  few 
published  studies.  The  results  of  these  determinations 
are  summarized  in  Table  I.  Included  are  results  from 
Chapter  III  of  this  thesis,  which  formed  part  of  a  recent 
publication  (62). 
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The  results  in  Table  I  and  the  corresponding  published 
comments  invite  several  conclusions  regarding  trimethyllead 
complexes.  With  the  exception  of  hydroxide  the  affinity 
of  these  ligands  for  trimethyllead  is  small  relative  to 
that  of  water.  Hence,  water  is  the  only  ligand  known  to 
form  two  to  one  ligand  to  metal  complexes  with  trimethyl¬ 
lead  in  aqueous  solution.  However,  a  second  trimethyllead 
cation  may  bind,  albeit  weakly,  to  the  coordinated  donor 
atom  of  a  relatively  strong  one-to-one  complex  such  as 
trimethyllead  hydroxide.  In  addition,  a  direct  correlation 
has  been  proposed  to  exist  between  trimethyllead  binding 
constants  and  ligand  basicity  (77) ;  in  the  absence  of 
other  factors,  strongly  basic  ligands  should  prove  to  be 
the  most  effective  complexing  agents  for  trimethyllead, 
a  Lewis  acid.  This  correlation  will  be  examined  in 
Chapter  IV  of  this  thesis. 

A  number  of  the  ligands  in  Table  I  are  present  in 
biological  fluids.  At  physiological  pH  levels,  there 
will  be  significant  fractions  of  trimethyllead  complexed 
by  each  of  these,  even  those  with  small  formation  con¬ 
stants  (63) .  The  presence  of  such  complexes  may  account 
in  part  for  the  unique  biological  properties  of  trimethyl¬ 


lead  (25)  . 
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C .  The  Present  Study 

Triorganolead  systems  have  been  studied  extensively 
in  recent  years  by  means  of  various  nmr  techniques  (44,49, 
54-57,59,62,63,67,69,77-83).  Several  of  these  studies 
have  demonstrated  that  the  nmr  chemical  shift  of  tri- 
methyllead  protons  provides  a  sensitive  probe  of  the 
extent  of  complex  formation  (62,63,77).  The  variation  of 
this  chemical  shift  with  pH  was  employed  in  this  thesis 
to  determine  the  equilibrium  formation  constants  of 
aqueous  trimethyllead  with  hydroxide  and  with  various 
organic  ligands.  Details  of  this  technique  will  be 
presented  in  the  appropriate  sections  of  this  thesis. 

Complexes  of  trimethyllead  with  ligands  containing 

sulfhydryl  groups  form  the  majority  of  the  species  studied. 

Several  factors  influenced  the  selection  of  ligands: 

sulfhydryl  groups  are  common  in  biological  systems,  they 

are  relatively  basic  when  deprotonated  and  such  ligands 

have  been  found  to  complex  very  strongly  with  methyl- 

mercury  (84) ,  the  central  metal  atom  of  which  is  iso- 

electronic  with  that  of  trimethyllead.  A  major  factor 

2+ 

in  choosing  sulfhydryl  ligands  is  that  Pb  is  strongly 
complexed  by  sulfhydryl  groups. 

This  work  marks  the  conclusion  of  a  program  of 
research  directed  towards  quantitatively  characterizing 
the  aqueous  coordination  chemistry  of  trimethyllead. 


CHAPTER  II 


EXPERIMENTAL 


A.  Chemicals 

Glycine  (AnalaR,  British  Drug  Houses  Chemicals  Ltd.), 
2-mercaptoethanol  (Eastman  Kodak  Co.),  N-acetyl-D , L- 
penicillamine  (Sigma  Chemical  Co.),  D , L-penicillamine 
(Aldrich  Chemical  Co.),  L-cysteine  (Nutritional  Biochemicals 
Corp. ) ,  reduced  glutathione  (Sigma  Chemical  Co.)  and 
t-butanol  (J.  T.  Baker  Chemical  Co.)  were  used  as  received. 
In  some  cases,  the  above  reagents  or  their  solutions  were 
standardized  as  described  below.  Trimethyllead  acetate 
(Ventron  Alfa  Products)  was  the  only  water  soluble  form 
of  trimethyllead  commercially  available.  To  avoid  inter¬ 
ference  due  to  acetate  complexation  of  trimethyllead 
(62,77),  the  trimethyllead  acetate  was  converted  to  a 
trimethyllead  perchlorate  stock  solution.  The  purification 
by  ion  exchange  and  the  standardization  by  pH  titration  of 
this  solution  is  described  below. 

All  other  chemicals  were  of  the  highest  grade  commer¬ 
cially  available  and  were  used  without  further  purification 
or  standardization.  Doubly  distilled  water  was  used 
exclusively  in  all  stages  of  solution  preparation  and  as 
a  final  rinse  in  the  cleaning  of  all  glassware. 
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B .  Preparation  and  Standardization  of  Trimethyllead 

Perchlorate  Solutions 

Ion  exchange  has  been  used  to  remove  acetate  from 
methylmercuric  hydroxide  solutions  (85)  and  a  similar 
method  was  developed  for  the  conversion  of  trimethyllead 
acetate  to  trimethyllead  perchlorate.  A  column  of  Dowex  2- 
X8  anion  exchange  resin  was  prepared  for  use  by  first 
washing  with  0.5  M  hydrochloric  acid  and  doubly  distilled 
water  to  remove  residual  impurities,  and  then  regenerating 
with  0.5  M  sodium  hydroxide.  A  freshly  prepared  solution 
of  approximately  0.25  M  trimethyllead  acetate  was  filtered 
under  vacuum  through  a  Nalgene  disposable  filter  unit  of 
0.2  ym  pore  size  to  remove  an  insoluble  fraction.  These 
filter  units  were  found  to  release  a  yellowish  impurity 
(possibly  a  membrane  binding  agent)  into  the  filtrate 
when  basic  solutions  were  filtered.  Therefore,  all  filter 
units  were  rinsed  several  times  before  use  by  drawing 
0.5  M  sodium  hydroxide,  followed  by  doubly  distilled  water, 
through  the  filter  until  the  filtrate  was  colorless. 

Proton  magnetic  resonance  spectra  of  the  resultant  solu¬ 
tions  showed  no  unexpected  resonances. 

The  filtrate  was  then  passed  several  times  through 
the  anion  exchange  column;  the  number  of  passes  required 
depended  on  the  volume  of  stock  solution  to  be  prepared 
and  on  the  capacity  of  the  column.  The  eluate  was  filtered 
through  fresh  Nalgene  filter  units  after  each  pass  to 
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remove  a  cloudy  white  precipitate  which  formed  soon  after 
elution.  This  precipitate  was  thought  to  be  trimethyllead 
hydroxide  (the  hydroxides  of  most  cations  are  less  soluble 
than  the  corresponding  acetates  and  in  addition,  the 
original  trimethyllead  acetate  solution  was  near  saturation 
before  the  exchange  of  ions  took  place  in  the  column) .  The 
e-luate  was  filtered  to  eliminate  the  possibility  of  co¬ 
precipitation  of  residual  trimethyllead  acetate.  Pre¬ 
cipitate  also  appeared  after  the  final  pass  through  the 
column,  but  dissolved  completely  below  pH  9.7  as  the  pH 
was  adjusted  with  perchloric  acid. 

The  column  was  regenerated  between  passes  with  0.5  M 
sodium  hydroxide.  The  proton  magnetic  resonance  spectrum, 
at  high  amplitude,  of  a  sample  was  taken  after  each  pass 
to  test  for  the  presence  of  residual  acetate.  These  pmr 
samples  were  first  acidified  with  concentrated  perchloric 
acid?  in  basic  solution  a  satellite  resonance  of  trimethyl¬ 
lead  partially  obscures  the  small  acetate  resonance. 

In  acidic  solution,  the  resonance  of  the  methyl  protons 
of  the  resultant  acetic  acid  is  located  approximately 
0.84  ppm  downfield  from  the  tert-butyl  resonance  of 
t-butanol  and  it  is  detectable  at  concentration  levels 
of  greater  than  approximately  0.001  M  with  the  nmr 
instrument  used  in  this  research.  When  the  acetic  acid 
resonance  fell  below  this  detection  limit  in  the  presence 
of  0.2  M  trimethyllead  perchlorate,  the  solution  was 


refiltered  through  a  fresh  Nalgene  unit,  neutralized  to 
between  pH  5  and  7 ,  and  stored  tightly  sealed  in  the 
absence  of  light  (21) . 

The  stock  solution  of  trimethyllead  perchlorate  was 
standardized  by  titration  with  sodium  hydroxide  using  a 
Mettler  DVll  Automatic  Titrator,  using  a  glass  electrode  to 
locate  the  endpoint.  To  ensure  complete  neutralization 
of  all  the  trimethyllead  hydroxide  present  in  solution, 
the  initial  pH  was  first  reduced  to  a  value  between  pH 
2  and  3  using  concentrated  perchloric  acid.  The  solution 
is  thus  a  mixture  of  a  strong  acid  (HCIO^)  and  a  weak  acid 
((CH^)3Pb+),  which  yields  a  titration  curve  with  two 
inflection  points.  Titration  of  the  added  perchloric  acid 
gives  a  sharp  inflection  point  at  pH  5  and  the  (CH^^Pb 
a  less  well. defined  inflection  point  at  approximately 
pH  10.4.  The  exact  position  of  the  trimethyllead  endpoint 
was  determined  using  a  second  derivative  plot  of  the 
titration  curve  (86) .  The  titration  with  base  converted 
the  trimethyllead  from  (CH3)3Pb+  to  (CH3)3PbOH.  Consequent¬ 
ly,  the  amount  of  sodium  hydroxide  added  between  the  first 
and  second  endpoints  provided  the  concentration  of  tri¬ 
methyllead  in  the  stock  solution. 

During  the  period  of  development  of  this  technique, 
the  existence  of  a  third  species  of  trimethyllead  in 
aqueous  solution  (i.e.,  the  dimer,  (  (CH3 ) 3Pb) 3OH  )  was 
suspected  (77) .  Verification  of  the  accuracy  of  the 
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po tent iome trie  technique  was  therefore  required.  This  was 
accomplished  by  determining  the  trimethyllead  concentration 
of  several  stock  solutions  with  atomic  absorption  spectro¬ 
scopy  (62,77),  and  with  EDTA  titrations  of  aliquots 

24- 

following  conversion  to  Pb  .  The  atomic  absorption 
results  tended  to  be  slightly  higher  than  the  trimethyl¬ 
lead  concentrations  obtained  by  potentiometry ,  but  the 
agreement  was  within  experimental  error,  supporting  the 
potentiometric  standardization.  Further  confirmation  of 
the  accuracy  of  the  potentiometric  method  was  obtained  with 
the  reduction-EDTA  titration  method,  which  also  gave 
results  identical  within  experimental  error  to  those  of 
potentiometry . 

The  procedure  used  in  the  reduction-EDTA  titration 

method  involved  the  digestion  of  aliquots  (typically  5.00 

ml)  of  the  stock  solution  with  15  ml  of  aqua  regia  (one 

part  concentrated  nitric  acid  plus  four  parts  concentrated 

hydrochloric  acid)  for  about  four  hours  to  ensure  complete 

4-  24- 

reduction  of  (CH^^Pb  to  Pb  .  The  sides  and  cover  of 
the  reaction  vessel  were  periodically  rinsed  down  with 
concentrated  nitric  acid.  The  solutions  were  then  evapo¬ 
rated  to  one  half  of  their  original  volumes.  Several 
volumes  of  nitric  acid  were  added  during  a  period  of  several 
hours  of  gentle  heating  in  order  to  evaporate  chloride 
(as  HCl(g))  which  causes  precipitation  of  PbCl2-  The 
resultant  acidified  Pb(NC>3)2  solutions  were  adjusted  to 
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approximately  neutral  pH  with  solid  sodium  hydroxide 
pellets.  The  pH  value  was  followed  with  a  combination 
glass  electrode.  This  solution  was  titrated  with  standar¬ 
dized  EDTA  using  hexamine  buffer  and  xylenol  orange  indi¬ 
cator  (87) .  The  EDTA  solution  was  standardized  by  titration 
with  a  ZnSO^  solution  of  accurately  known  concentration  in 
t-he  presence  of  xylenol  orange  indicator  (87)  . 

It  was  necessary  to  modify  slightly  Vogel's  (87) 

2+ 

procedure  for  the  titration  of  Pb  .  The  quantity  of 
hexamine  specified  did  not  always  provide  sufficient  buffer 
capacity  to  prevent  a  decrease  in  pH  such  that  suboptimal 
levels  were  reached  during  the  course  of  the  titration. 

This  drop  in  pH  is  due  to  the  release  of  protons  from  di¬ 
sodium  EDTA  during  complexation .  Continuous  monitoring 
of  pH  with  a  small  combination  glass  electrode,  together 
with  further  addition  of  solid  sodium  hydroxide  as  required 
to  maintain  a  pH  of  about  6,  ensured  consistently  good 
indicator  endpoints. 

The  three  independent  procedures  yield  identical  re¬ 
sults  within  experimental  error,  supporting  the  validity 
of  these  methods.  The  presence  of  a  small  amount  of  the 
dimeric  species,  (  (CH^ )  ^Pb)  20H+ ,  was  not  f°un<3  to  interfere 
in  the  pH  titrations  at  the  concentrations  used  here. 

(It  is  now  known  (and  detailed  in  Chapter  III)  that  this 
species  is  significant  only  in  the  neutral  pH  regions, 
and  readily  reacts  with  a  further  hydroxide  ion.  At  the 
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end  point,  all  trimethyllead  in  solution  is  present  as 
trimethyllead  hydroxide;  thus,  a  one  to  one  stoichiometry 
will  be  preserved  in  titrations  of  (CH^J^Pb*  with  hydroxide). 
Since  the  potentiometric  method  is  quite  convenient  and 
yields  the  most  precise  results,  it  was  the  method  of 
choice  for  the  standardization  of  succeeding  preparations 
of  trimethyllead. 

C.  Preparation  and  Standardization  of  Ligand  Solutions 

All  solutions  in  this  and  subsequent  sections  were 
prepared  with  doubly  distilled  water  which  had  been  freshly 
boiled  and  cooled  under  a  stream  of  argon  (to  prevent  the 
reabsorption  of  oxygen  and  carbon  dioxide) .  All  pipets 
were  calibrated  by  a  standard  method  (89) . 

Liquid . 2-mercaptoethanol  (nominally  98%  pure)  was 
the  source  of  the  mercaptoethanol  ligand.  A  5.00  ml  aliquot 
of  the  liquid  was  diluted  to  100  ml  with  water.  A  2.00  ml 
aliquot  of  the  dilute  solution  was  added  to  140  ml  of  water, 
to  which  sufficient  crystalline  sodium  perchlorate  had 
been  added  to  give  a  final  ionic  strength  (y)  of  0.3  M.  The 
concentration  of  ligand  in  the  final  solution  was  approxi¬ 
mately  0.012  M.  Concentrated  perchloric  acid  was  added 
dropwise  until  the  pH  of  the  solution  was  between  2.6  and 
3.0.  The  added  acid  contributed  only  about  0.01  M  to  the 
final  ionic  strength.  The  resultant  solution  was  titrated 
with  standard  0.12  M  sodium  hydroxide  using  a  Mettler 


DVll  Autotitrator .  Titrant  was  added  in  small  increments 
in  the  regions  of  the  endpoints  and  near  the  half-neutrali 
zation  point  of  2-mercaptoethanol  to  define  these  regions 
as  well  as  possible.  The  titration  curve  was  followed 
with  a  glass  electrode.  As  with  the  standardization 
of  trimethyllead ,  the  solution  was  a  mixture  of  a  strong 
acid  (perchloric  acid)  and  a  weak  acid  (2-mercaptoethanol) 
thus  giving  two  inflection  points  in  the  curve.  The 
sharply  defined  inflection  point  at  pH  5  was  due  to  HCIO^ 
while  the  less  distinct  break  at  pH  10.6  was  due  to  the 
ligand.  The  positions  of  the  endpoints  were  determined 
using  a  second  derivative  plot  of  the  titration  curve  (86) 
The  amount  of  sodium  hydroxide  added  between  the  two 
endpoints  provided  the  concentration  of  2-mercaptoethanol 
in  the  dilute  solution,  and  hence,  in  the  concentrated 
liquid.  Ten  replicates  gave  a  median  purity  of  97.2  % 

±  0.6%. 

A  proton  magnetic  resonance  integration  experiment 
was  performed  to  determine  the  species  present  in  the 
concentrated  liquid  2-mercaptoethanol.  Six  replicates  of 
the  integration  curves  of  the  neat  liquid  2-mercapto¬ 
ethanol  spectrum  were  obtained  on  a  Varian  HA-100  spectro¬ 
meter.  Ratios  of  the  heights  of  the  integration  waves 
were  compared,  and  the  results  obtained  were  consistent 
with  a  solution  composed  of  97%  HOCH2CH2SH,  1.5% 
HOCH2CH2SSCH2CH2OH  and  1.5%  I^O.  Neither  impurity  inter- 
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feres  with  the  nmr  determination  of  the  formation  constants 
for  the  complexation  of  (CH^J^Pb"1"  by  the  sulfhydryl  group 
of  2-mercaptoethanol ,  provided  that  standardized  solutions 
are  used.  The  nmr  integration  experiment  also  confirms 
the  accuracy  of  the  pH  titration  standardization.  The 
determination  of  the  acid  dissociation  constant  of  the 
ligand  thiol  group  from  the  titration  data  is  described 
in  section  I  of  this  chapter. 

N-acetyl-D ,L-penicillamine  was  standardized  using  a 
method  similar  in  most  respects  to  that  used  for  2-mercapto- 
ethanol.  The  significant  differences  were  that  the  ligand 
was  added  directly  to  the  titration  mixture  as  a  crystalline 
solid  and  no  perchloric  acid  was  required.  The  ligand 
carboxylic  acid  group  provided  a  good  endpoint  at  pH  7, 
while  the  sulfhydryl  group  endpoint  was  of  no  use,  being 
poorly  defined  in  the  high  pH  region.  To  obtain  the  desired 
information,  it  was  assumed  that  the  concentration  of  car¬ 
boxylic  acid  groups  equalled  that  of  the  sulfhydryl  groups. 
Hence,  the  volume  of  standard  sodium  hydroxide  required 
to  titrate  to  the  first  endpoint  was  used  to  calculate  the 
purity  (in  terms  of  the  sulfhydryl  group)  of  the  crystalline 
solid.  Four  replicates  gave  a  median  value  of  99.1  (±  0.4) 
These  results  were  corroborated  by  several  determinations 
of  sulfhydryl  purity  using  the  standardization  procedure 
described  by  Benesch  and  Benesch  (90) .  In  this  procedure, 
the  reaction  of  the  reduced  ligand  thiol  groups  with 
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iodoacetamide  stoichiometrically  releases  into  solution 
the  protons  of  all  such  groups.  The  concentration  of  the 
protons  released  is  determined  by  appropriate  titrimetric 
measurements  (90).  This  procedure  gave  a  median  result  of 
98.7  (±1.0)  %,  which  confirms  the  accuracy  of  the  pH 
titration  technique  within  experimental  error. 

The  purity  of  the  reduced  sulfhydryl  group  of  the 
glutathione  employed  in  this  study  was  determined  by 
titration  with  coulometrically  generated  iodine  in  acetate 
buffer  using  biamperometric  endpoint  detection  (91) .  Five 
replicates  gave  a  purity  with  respect  to  the  -SH  group 
of  97.3  (±  0.5)  %  (93).  The  purities  of  the  -SH  groups 
of  L-cysteine  and  D,L-penicillamine  were  determined  by 
titration  with  coulometrically  generated  bromine  in  1  M 
hydrochloric  acid  using  biamperometric  endpoint  detection 
(92).  Five  replicates  gave  a  purity  for  cysteine  of  95.0% 
±  0.5%  and  for  penicillamine  of  99.6%  ±  0.4%  (93). 

The  sodium  hydroxide  solutions  used  in  all  phases  of 
this  work  were  prepared  under  an  atmosphere  of  argon  using 
the  carefully  degassed  water  described  at  the  beginning  of 
this  section.  Standardization  of  these  solutions  was 
accomplished  by  titration  against  known  amounts  of  primary 
standard  potassium  hydrogen  phthalate  (88) .  The  standard 
sodium  hydroxide  stored  in  the  autotitrator  reservoir  was 
protected  from  atmospheric  carbon  dioxide  by  means  of  an 
Ascar ite-containing  tube  on  the  only  air  inlet.  This 
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stock  solution  was  restandardized  carefully  before  each 
use,  even  though  its  molarity  was  found  to  be  stable  for 
several  weeks. 

D.  pH  Measurements 

pH  Measurements  were  carried  out  at  25  ±  1  °C  with 
an  Orion  Model  701  digital  pH  meter  equipped  with  a  standard 
full  pH  range  glass  electrode  and  a  silver/silver  chloride 
reference  electrode  in  which  saturated  sodium  sulfate  was 
used  as  the  electrolyte  solution.  This  reference  electrode 
was  used  rather  than  the  standard  Ag/AgCl  reference 
electrode  with  potassium  chloride  electrolyte  to  avoid 
problems  resulting  from  the  precipitation  of  potassium 
perchlorate  in  the  reference  electrode  fiber  junction. 

Fisher  Certified  Buffer  solutions  having  pH  values  of  4.01, 
7.00  and  10.00  at  25°C  (potassium  hydrogen  phthalate, 
potassium  phosphate  and  potassium  borate/potassium  carbonate, 
respectively)  were  used  to  standardize  the  pH  meter. 

E .  Proton  Nuclear  Magnetic  Resonance  Measurements 
Proton  nmr  spectra  were  obtained  on  a  Varian  A-60-D 

high  resolution  spectrometer  at  a  probe  temperature  of 
25  ±  1 °C .  The  temperature  of  the  probe  was  determined 
by  measuring  the  potential  of  a  copper  vs.  constantan 
thermocouple  inserted  in  the  probe. 

Spectra  were  recorded  at  sweep  rates  of  0.2  Hz/sec 


27 


for  the  chemical  shift  measurements.  Reported  data  are 
the  median  of  three  scans.  Chemical  shifts  were  measured 
relative  to  the  tert-butyl  resonance  of  t-butanol  (TBA) 
or  the  central  resonance  of  tetramethylammonium  (TMA)  cation, 
but  are  reported  relative  to  the  methyl  resonance  of 
sodium  2 , 2-dimethyl-2-silapentane-5-sulf onic  acid  (DSS) . 
Positive  shifts  correspond  to  resonances  of  protons  less 
shielded  than  the  methyl  protons  of  DSS.  The  tert-butyl 
resonance  of  TBA  is  1.243  ppm  downfield  from  the  methyl 
resonance  of  DSS;  the  central  resonance  of  TMA  is  3.175 
ppm  downfield  from  DSS. 

The  sweep  widths  of  the  nmr  recorder  were  calibrated 
by  a  sideband  procedure.  A  Hewlett-Packard  Model  200  AB 
audio  oscillator  was  used  to  provide  sideband  signals 
from  tetramethylsilane  (12%  in  chloroform) .  The  frequency 
of  the  sideband  was  counted  with  a  Hewlett-Packard  Model 
5307A  High  Resolution  Counter.  For  recorder  calibration, 
the  shift  of  the  sideband  from  the  main  signal  (as  indicated 
by  the  calibration  lines  of  the  recorder  paper)  was  matched 
to  the  counted  frequency  generated  by  the  audio  oscillator. 

F.  Carbon-13  Nuclear  Magnetic  Resonance  Measurements 

Carbon-13  nmr  spectra  were  obtained  using  a  Varian 
HA-100  spectrometer  operating  at  a  frequency  of  25.1  MHz 
and  equipped  with  a  Digilab  FTS/NMR-3  Data  System.  The 
pulsed  Fourier  transform  mode  was  used  with  proton  de- 
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coupling.  For  each  free  induction  decay  signal,  8  K  data 

points  were  collected  in  the  computer.  The  spectra  were 

time  averaged  using  500  accumulations  in  all  cases  except 

one,  where  4  K  accumulations  were  required  to  measure  the 
2  07_,  13n 

Pb  -  C  coupling  constant  at  pH  5.0.  These  accumu¬ 
lations  gave  adequate  signal-to-noise  ratios  in  each 
case.  The  deuterium  resonance  from  in  a  coaxial 

D  D 

capillary  served  as  the  lock  signal.  The  frequency  range 
of  the  transformed  spectra  was  5000  Hz.  Chemical  shifts 
are  reported  in  ppm  relative  to  the  carbon-13  resonance 
of  the  1,4-dioxane  added  as  an  internal  reference.  The 
resonance  of  1,4-dioxane  is  67.73  ppm  downfield  from  the 
tetramethylsilane  carbon-13  resonance.  Positive  chemical 
shifts  correspond  to  less  shielding  than  in  1,4-dioxane. 

The  chemical  shift  measurements  are  considered  accurate  to 
within  0.05  ppm.  For  all  measurements  the  sample  tempera¬ 
ture  was  25  -  1°C  as  measured  with  a  standard  thermocouple. 

G.  Solutions  for  Nuclear  Magnetic  Resonance  Measurements 

Solutions  used  in  the  nuclear  magnetic  resonance 
measurements  were  prepared  in  carefully  degassed  doubly 
distilled  water  from  buret-delivered  aliquots  of  stock  tri- 
methyllead  perchlorate.  The  ligand  was  weighed  into  the 
mixture  as  a  solid  or  an  aliquot  of  ligand  solution  was  added 
to  the  mixture  by  buret  or  pipet.  The  concentration  of 
trimethyl lead  in  the  solutions  from  which  nmr  measurements 
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were  determined  was  0.005  M  to  0.075  M.  Ratios  of  ligand 
to  trimethyllead  varied  between  0.5  and  10  to  1,  but 
typically  were  either  2  or  4  to  1.  Tert-butanol  was 
added  as  an  internal  chemical  shift  reference  for  the  tri¬ 
methyllead  resonance,  at  a  concentration  about  two-thirds 
that  of  trimethyllead.  Tetramethylammonium  nitrate  was 
used  as  the  internal  standard  instead  of  TBA  only  for  those 
ligand  shift  measurements  where  TMA  allowed  the  use  of  a 
narrower  nmr  sweep  width.  The  concentration  of  TMA  was 
typically  less  than  0.01  M.  In  all  cases,  the  ionic  strength 
was  maintained  at  approximately  0.3  M  with  sodium  perchlorate. 

In  the  experiments  which  measured  chemical  shift  as  a 
function  of  pH,  the  ionic  strength  of  the  solution  was 
first  adjusted  to  0.3  M  and  then  samples  were  withdrawn 
as  the  pH  was  varied.  Variations  in  ionic  strength  caused 
by  change  in  the  chemical  species  in  solution  are  small 
(typically  t  0.02  M)  relative  to  the  background  ionic 
strength.  The  solutions  were  initially  made  acidic  to 
about  pH  2.5,  using  a  measured  amount  of  perchloric  acid. 

Then  sodium  hydroxide  was  added  and  samples  were  withdrawn 
at  appropriate  pH  intervals  until  approximately  pH  12.5 
was  reached.  Usually  twenty-one  to  twenty-five  0.4  ml 
samples  were  required  to  cover  this  pH  range.  Concentrated 
perchloric  acid  (70%)  and  deoxygenated  carbonate-free 
sodium  hydroxide  solutions  were  used  to  adjust  the  pH. 

Since  compounds  containing  thiol  groups  are  subject  to 


air  oxidation,  and  since  high  pH  samples  absorb  atmos¬ 
pheric  carbon  dioxide  (causing  downward  pH  drift) ,  all 
possible  stages  of  the  solution  preparation  and  sampling 
were  performed  under  a  stream  of  argon.  Nmr  spectra  were 
obtained  immediately  after  sample  preparation  or  after 
overnight  storage  of  samples  under  argon  at  1°C. 

In  the  mole  ratio  experiment  of  trimethyllead  vs. 
2-mercaptoethanol ,  each  sample  solution  was  made  separately 
with  varying  ligand  to  trimethyllead  ratios,  but  at 
constant  pH  and  ionic  strength.  The  ratio  of  ligand  to 
trimethyllead  in  these  solutions  ranged  from  0.0  to  6.1. 

The  ionic  strength  of  all  samples  was  held  at  approximately 
0.3  M  with  sodium  perchlorate.  The  pH  in  all  cases  was 
adjusted  to  8.0;  very  little  base  was  added. 

The  solutions  for  carbon-13  nuclear  magnetic  resonance 
spectroscopy  were  prepared  in  a  manner  similar  to  that  for 
the  proton  nmr  experiments,  except  that  dioxane  was  used 
as  the  internal  standard.  The  concentration  of  dioxane 
was  about  0.04  M,  the  trimethyllead  was  0.0947  M,  the 
glutathione  was  0.146  M  and  the  perchlorate  was  approxi¬ 
mately  0.17  M  (added  as  70%  HCIO^)  giving  a  y  of  0.4  M. 

H.  Determination  of  Formation  Constants  by  NMR  Spectro¬ 
scopy 

The  formation  constants  of  the  trimethyllead  complexes 


were  calculated  from  the  chemical  shift  of  the  methyl 


< 
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protons  of  trimethyllead  in  solutions  of  known  pH  and  known 
trimethyllead  and  ligand  concentrations. 

All  the  systems  studied  in  this  thesis  were  found 
to  be  labile;  that  is,  the  trimethyllead  moiety  is  ex¬ 
changing  between  its  various  forms  (aquated,  hydroxy- 
complexed  and  ligand-complexed)  at  a  rate  which  is  fast  on 
the  nmr  time  scale.  Thus,  a  single  sharp  resonance  is 
observed  for  the  methyl  protons  of  trimethyllead.  The 
observed  chemical  shift  of  this  resonance  is  the  weighted 
average  of  the  chemical  shifts  of  the  various  trimethyllead 
species.  This  is  expressed  by  the  following  relation: 

«OBS  =  £  Vi 

where  <5qBS  represents  the  observed,  exchange  averaged 

chemical  shift  of  trimethyllead,  P.  is  the  fraction  of  the 
th 

i  trimethyllead  species  in  solution  and  6^  is  the  chemical 
shift  of  that  species. 

The  formation  constants  and  the  chemical  shifts  of 
the  complexes  were  obtained  from  pH  vs  chemical  shift 
experiments  by  curve-fitting  the  data  to  Equation  1.  The 
curve-fitting  program  used  was  generously  supplied  by 
Drs.  J.  L.  Dye  and  V.  A.  Nicely  of  Michigan  State  Univer¬ 
sity.  This  program,  KINET,  has  been  described  in  the 
literature  (94).  A  user-supplied  subroutine  is  required 
to  describe  the  model  to  be  fitted  to  the  data.  To 
use  Equation  1,  it  is  first  necessary  to  solve  for  the 
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various  fractions,  P.^.  The  fractions  are  calculated  from 
the  equations  for  the  formation  constants  and  the  mass 
balance  equations  for  trimethyllead  and  ligand.  Manipula¬ 
tion  of  the  i  equations  to  eliminate  all  fractions  but 

one,  P  ,  gives  a  polynomial  in  P  .  The  coefficients  of 
x  x 

P  in  this  polynomial  are  in  terms  of  the  formation  con- 

A 

stants,  the  concentrations  of  hydrogen  ion  or  hydroxide 
ion,  and  the  total  concentrations  of  trimethyllead  and  of 
ligand.  Solution  of  this  polynomial  for  P  was  achieved 

X 

iteratively  using  the  Library  Subroutine  DRTMI ,  which  is 
part  of  the  IBM  Scientific  Subroutine  Package  (SSP  Library) . 
With  this  value  for  P  at  a  particular  pH,  the  other  frac- 

X 

tions  were  calculated  using  the  appropriate  formation 
constant  equations.  The  numerical  values  for  these 
fractions  were  used  to  calculate  a  chemical  shift  with 
Equation  1.  This  entire  process  was  repeated  for  each  data 
point.  The  sum  of  the  squares  of  the  differences  between 
each  predicted  chemical  shift  and  the  corresponding  observed 
chemical  shift  was  minimized  by  the  program  KINET.  It 
accomplished  this  by  systematically  varying  the  initial 
estimates  for  the  unknown  formation  constants  and  chemical 
shifts  until  the  predicted  and  observed  curves  converged. 

The  values  of  the  unknown  parameters  at  convergence  were 
considered  accurate  only  if  all  the  error  indicators 
calculated  by  KINET  were  small.  In  order  to  determine 
accurate  formation  constants  by  the  above  method,  it  was 
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necessary  to  have  in  solution  a  significant  fractional 
concentration  of  each  form  of  trimethyllead  in  some  portion 
of  the  observed  pH  range.  Those  complexes  which  did  not 
meet  this  criterion,  or  whose  chemical  shift  was  similar 
to  that  of  a  species  in  plentiful  supply,  received  from 
KINET  large  error  estimates  for  their  best-fit  formation 
constants  and  chemical  shifts. 

Examples  of  the  above  procedure  are  described  in  detail 
in  Chapters  III  and  IV. 

I .  Determination  of  Acid  Dissociation  Constants  of 
Selected  Ligands 

1  „  Determination  by  NMR  Measurements 

It  is  well  known  that  the  shifts  of  carbon-bonded 

protons  close  to  acidic  functional  groups  are  sensitive 

to  the  deprotonation  of  that  functional  group  (95) :  the 

method  outlined  in  the  previous  section  is  applicable  to 

the  determination  of  K  for  organic  acids  (96) .  The  acid 

a 

dissociation  constants  of  selected  ligands  were  calculated 
from  the  pmr  chemical  shift  of  the  ligand  resonances  in 
solutions  of  known  pH  and  ligand  concentrations.  The 
method  is  similar  to  that  described  above  for  the  deter¬ 
mination  of  formation  constants,  with  only  minor  differences. 
Tetramethylammonium  nitrate  (0.002  M)  was  used  as  the 
internal  standard.  Somewhat  higher  concentrations  of 
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ligand  (approximately  0.15  M)  were  required  due  to  the 
lower  intensity  of  the  ligand  methylene  resonances  compared 
to  the  degenerate  methyl  resonances  of  trimethyllead . 
Therefore,  the  ionic  strength  was  somewhat  more  difficult 
to  maintain  at  the  desired  value  of  0.3  M.  The  presence 
of  multiple  ligand  resonances  with  chemical  shifts  sensitive 
to  deprotonation  allowed  separate  determinations  of  the 
acid  dissociation  constant  with  each  resonance.  These 
results  were  averaged  to  give  the  value  used  in  this  work. 

Further  details  of  each  acid  dissociation  constant 
determination  are  described  in  Chapter  IV. 

2 .  Determination  by  pH  Titration 

The  acid  dissociation  constants  of  two  of  the  di¬ 
basic  ligands  were  obtained  by  potentiometric  titration 
using  a  glass  electrode.  The  titrations  were  carried  out 
under  argon  using  degassed  sodium  hydroxide  as  detailed 
in  the  previous  sections  in  this  chapter.  The  ionic 
strength  was  controlled  at  0.3  M  with  crystalline  sodium 
perchlorate . 

Both  ligands  titrated  were  dibasic  acids,  but  the 

estimated  pK  values  were  quite  different,  thus  warranting 
^  a 

different  experimental  procedures.  2-Mercaptoethanol  was 
thought  to  have  a  pK1  of  approximately  9  and  a  pK2  of 
above  14.  It  may  be  argued  that  such  a  high  value  for 
pj<2  effectively  renders  2— mercaptoethanol  a  monobasic 


. 
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acid  in  water;  hence,  pl^  was  ignored.  The  procedure 
was  therefore  as  follows.  Initially,  the  pH  was  reduced 
to  below  3  using  concentrated  perchloric  acid.  Titration 
of  the  perchloric  acid  yields  a  sharp  inflection  point  at 
pH  5.5  and  titration  of  the  sulfhydryl  group  of  the  ligand 
gives  a  less  well-defined  inflection  point  at  about  pH 
10.6.  The  acid  dissociation  constant,  which  is  defined 
by  Equations  2  and  3,  was  determined 

H_L  t  H+  +  HL~  (2) 

[H+] [HL_] 

K,  =  -  (3) 

[h2l] 

from  the  pH  data  between  the  two  inflection  points. 

N-acetyl-D ,L-penicillamine  was  thought  to  have  a  pK^ 
of  about  3  and  a  pl^  of  approximately  10.  Therefore,  no 
perchloric  acid  was  added:  a  weighed  amount  of  dry  crystal¬ 
line  ligand  was  dissolved  and  titrated  as  detailed  in 
section  C  above.  and  ,  as  defined  by  Equations  2  to 

5  inclusive, 

HL“  t  H+  +  L2"  (4) 

[H+] [L2"] 

K9  =  - - -  (5) 

[HL  ] 

were  determined  from  the  pH  titration  data  between  the  start 
and  the  first  inflection  point,  and  between  the  two  inflec- 
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tion  points,  respectively.  In  all  cases,  the  exact  posi¬ 
tions  of  the  endpoints  were  determined  by  second  derivative 
plots  (86) .  All  pH  titration  experiments  were  designed 
to  minimize  dilution  errors  and  activity  corrections, 
using  the  methods  described  by  Albert  and  Serjeant  (97) . 


CHAPTER  III 


TRIMETHYLLEAD  SPECIES  AND  EQUILIBRIA  IN  AQUEOUS  SOLUTION 


At  the  initiation  of  this  research,  no  detailed  infor¬ 
mation  was  available  regarding  the  stability  of  trimethyl- 
lead  in  aqueous  solution  or  the  nature  of  the  aqueous 
acid-base  equilibria  of  trimethyllead .  In  section  A  of 
this  chapter,  the  results  of  a  long-term  study  of  the 
stability  of  aqueous  trimethyllead  at  various  pH  levels 
are  reported.  In  section  B,  nuclear  magnetic  resonance 
results  are  described  for  the  aqueous  acid-base  solution 
chemistry  of  trimethyllead.  The  aqueous  species  are  identi¬ 
fied  and  constants  for  the  equilibria  between  the  various 
species  are  determined.  The  research  described  in  section 
B  was  performed  in  collaboration  with  Dr.  T.  L.  Sayer, 

E.  K.  Millar  and  Dr.  C.  A..  Evans  (62). 

A.  The  Stability  of  Trimethyllead  in  Aqueous  Solution 

1 .  Introduction 

In  aqueous  solution,  the  methyllead  compounds 
(CH3)3PbX  and  (CH3)2PbX2  redistribute  according  to  the 
following  reaction  scheme  (18) : 


4 (CH3) 3PbX  t 

2(CH3)2PbX2  +  2(CH3)4Pb 

(6) 

2  (CH3)  2PbX2 

(CH3)3PbX  +  CH3X  +  PbX2 

(7) 

3 (CH3) 3PbX  + 

2(CH3)4Pb  +  CH3X  +  PbX2 

(8) 
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Reaction  7  is  much  faster  than  6  and  the  rate  of  both 
reactions  depends  on  the  nature  of  the  anion  X  (98,99). 

The  perchlorate  anion  used  extensively  in  the  present  study 
was  shown  to  produce  a  much  lower  rate  constant  for  reaction 
7  than  did  many  common  anions  (98) ,  especially  under  the 
conditions  employed  in  this  work.  The  overall  stoichio¬ 
metry  of  the  disproportionation  of  trimethyllead  can  be 
described  by  reaction  8. 

The  comparatively  low  rate  constant  of  reaction  6 
seems  to  indicate  that  trimethyllead  (and  especially  tri¬ 
methyllead  perchlorate)  is  relatively  stable  at  room  tem¬ 
perature  (18) .  Nonetheless,  the  occurrence  of  reactions  6 
through  8  as  well  as  other  observations  (52)  cast  some  doubt 
on  the  exact  composition  and  practical  stability  of  the 
trimethyllead  perchlorate  stock  solutions  used  in  this 
work.  Therefore,  a  study  to  determine  the  long  term 
stability  of  the  lead-carbon  bonds  towards  hydrolysis  at 
several  pH  levels  was  completed. 

2 .  Results  and  Discussion 

Solutions  of  0.185  M  trimethyllead  at  pH  levels  of 
0.9,  2.0,  7.1  and  11.0  were  tested  periodically  for  the 
presence  of  Pb(II)  ions  using  chemical  spot  tests.  Sodium 
rhodizonate  was  used  as  detailed  by  Feigl  and  Anger  (100) ; 

this  method  was  sensitive  to  inorganic  lead  at  concentrations 

-4 

greater  than  approximately  10  M. 


Since  the  indicator 
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was  inactive  at  high  pH,  the  alkaline  sample  spots  were 

neutralized  before  testing  using  the  tartrate  buffer 

described  in  the  reference  above.  Quantitative  results 

were  obtained  by  comparing  the  test  results  with  those  of 

a  fresh  series  of  standard  lead (II)  nitrate  solutions, 

-1  -5 

which  ranged  in  concentration  from  10  M  to  10  M. 

After  28  months  of  monitoring,  the  results  indicated  that 

-4 

the  concentration  of  inorganic  lead  was  less  than  10  M 

in  solutions  at  neutral  pH.  This  corresponded  to  less  than 

one  part  per  thousand  decomposition  of  trimethyllead 

perchlorate  to  inorganic  lead.  Lead (II)  ions  were  found 

-2 

to  be  present  at  levels  of  approximately  10  M  in  pH 

-2  -3 

0.9  solutions,  between  10  M  and  10  M  in  pH  2 . 0  solu- 

-3  -4  . 

tions  and  between  10  and  10  M  in  pH  11.0  solutions. 

These  results  were  corroborated  by  the  comparison  of 
the  integrated  proton  nmr  intensities  of  the  above  samples 
during  the  initial  five  month  experimental  period.  No 
new  resonances  appeared  as  time  progressed.  In  addition, 
the  intensity  measurements  of  the  methyl  resonance  of 
trimethyllead  relative  to  those  of  the  internal  standard 
tert-butanol  were  consistent  with  the  spot  test  results 
within  experimental  error.  Neutral  pH  was  thus  found  to 
be  the  most  favorable  for  prolonged  stability  of  tri¬ 
methyllead  perchlorate  in  aqueous  solution. 

Taken  together,  these  results  established  that  aqueous 
solutions  of  trimethyllead  perchlorate  do  not  decompose 
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significantly,  even  over  long  periods  of  time,  under  the 
conditions  employed  in  the  above  experiments.  Hence,  these 
conditions  were  met  in  the  storage  of  all  trimethyllead 
stock  solutions:  the  maintenance  of  darkness  under  an 
inert  atmosphere  (14)  at  neutral  pH,  the  use  of  the  per¬ 
chlorate  anion  only  and  the  use  of  concentrations  of  tri¬ 
methyllead  below  the  observed  maximum  of  about  0.2  M.  For 
shorter  periods  of  time,  certainly  on  the  time  scale  of 
the  experiments  described  later  in  this  thesis,  trimethyl¬ 
lead  in  aqueous  solution  appeared  to  be  stable  to  hydrolysis 
of  the  lead-carbon  bond  even  at  extremes  of  pH. 

B.  The  Acid-Base  Chemistry  of  Trimethyllead 


1 .  Introduction 

The  proton  nuclear  magnetic  resonance  spectrum  of 

trimethyllead  in  aqueous  solution  consists  of  a  methyl 

singlet  flanked  symmetrically  by  two  less  intense  satellite 

peaks.  Figure  1  shows  a  typical  spectrum  at  approximately 

neutral  pH.  The  singlet  is  due  to  the  protons  of  the  methyl 

groups  bound  to  those  lead  atoms  having  a  nuclear  spin  of 

zero,  while  the  satellites  are  due  to  the  protons  of  the 

2  07 

methyl  groups  bonded  to  Pb,  the  isotope  having  a 
nuclear  spin  of  one— half  (22.6%  natural  abundance) . 

The  observed  chemical  shift,  ^qbS7  9^ven  ky  the 
position  of  the  central  resonance,  measured  from  the 
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resonance  of  a  chemical  shift  standard.  The  lead-proton 
spin-spin  coupling  constant,  J207pb-lH'  is  9iven  hY  the 
separation  of  the  satellite  lines.  Both  the  chemical 
shift  of  the  methyl  singlet  and  the  coupling  constant  are 
dependent  upon  the  nature  of  the  ligands  coordinated  to 
the  trimethyllead  in  solution.  In  aqueous  solutions  of 
trimethyllead  containing  no  coordinating  ligands  other 
than  water  and  hydroxide  ion,  the  chemical  shift  and  the 
coupling  constant  are  both  pH  dependent.  The  chemical  shift 
of  the  central  resonance  of  trimethyllead  perchlorate  in 
aqueous  solution  at  a  concentration  of  0.0100  M  is  shown  as 
a  function  of  pH  in  Figure  2.  (The  ionic  strength  was 
maintained  at  0.3  M  with  NaClO^.)  The  coupling  constant 
of  trimethyllead  varies  from  77  Hz  to  79  Hz  over  a  similar 
pH  range  (62,77).  The  pH  dependence  is  due  to  the  formation 
of  trimethyllead  hydroxide  complexes  as  the  concentration 
of  hydroxide  ion  is  increased.  No  evidence  exists  for 
the  coordination  of  more  than  one  hydroxide  ion  per  tri¬ 
methyllead  cation  up  to  at  least  pH  12.  Indeed,  the 
available  data  indicate  that  essentially  all  the  trimethyl¬ 
lead  is  present  as  (CH^^PbOH  at  high  pH. 

2  .  Results  and  Discussion 

The  aqueous  species  of  trimethyllead  and  the  equilibrium 
constants  for  these  species  were  determined  from  chemical 
shift  vs.  pH  data  of  the  type  displayed  in  Figure  2.  The 


I  I 
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Figure  2.  pH  dependence  of  the  chemical  shift  of  the  methyl 
protons  of  trimethyllead  in  a  0.0100  M  aqueous  solution  at 
0.3  M  ionic  strength  (maintained  with  NaClO^). 
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experimental  solutions  possessed  total  trimethyllead 
concentrations  ranging  from  5  x  10  M  to  2  x  10  M. 

At  each  pH  of  a  solution,  the  accepted  chemical  shift 
value  was  commonly  the  median  of  three  measurements. 

Lead-proton  spin-spin  coupling  constant  data  were  not 
used  to  determine  the  trimethyllead  species  and  equili¬ 
brium  constants.  The  small  relative  change  in  the  coupling 
constant  over  the  entire  range  of  pH  provided  data  of 
poor  sensitivity  towards  ligand  complexation :  accurate 
results  could  not  be  obtained  from  such  information.  These 
data  contrasted  with  the  data  obtained  from  the  analogous 
methylmercury  system  (85).  Both  the  mercury (199) -proton 
spin-spin  coupling  constants  and  the  chemical  shift  data 
could  be  used  to  obtain  precise  results,  because  both 
were  quite  sensitive  to  ligand  complexation. 

(a)  Determination  of  Species  Present  in  Aqueous 
Trimethyllead  Solutions.  Initially,  it  was  assumed  that 
a  simple  acid-base  equilibrium  existed  between  (CH3)3Pb+ 
and  (CH^J^PbOH,  and  that  no  other  species  of  trimethyllead 
was  present.  This  system  is  represented  by  the  following 
equations : 

(CH3)3Pb+  +  OH”  t  (CH3)3PbOH  (9) 

[ (CH-J .PbOH] 

K1  =  - — 

[ (CH3) 3Pb  ]aQH_ 


(10) 


where  aQH_  is  the  hydroxide  ion  activity  and  K1  is  the  for¬ 
mation  constant  of  (CH^J^PbOH.  The  limiting  chemical  shift 
at  a  pH  of  less  than  2  was  taken  to  be  the  shift  of  pure 
(CH^J^Pb  and  that  at  a  pH  of  greater  than  12  was  taken  to 
be  the  shift  of  pure  (CH^J^PbOH.  in  dilute  solutions  of 
trimethyllead ,  this  model  agreed  well  with  the  observed 
chemical  shift  data.  However,  the  displacement  of  the 
chemical  shift  curve  for  solutions  of  higher  trimethyllead 
concentration,  especially  in  the  neutral  and  moderately 
basic  pH  regions,  indicated  the  presence  of  a  third  species 
of  intermediate  chemical  shift.  In  addition,  this  con¬ 
centration  dependence  suggested  that  the  equilibria  between 
this  third  species  and  the  two  above-mentioned  species 
were  asymmetric;  the  position  of  an  asymmetric  equilibrium 
would  indeed  favor  the  simpler  species  (CH^J^Pb*  and 
(CH^ ) ^PbOH  at  lower  total  trimethyllead  concentrations. 

The  aqueous  methylmercury  system,  which  has  been 
described  in  detail  (85,101),  was  thought  to  be  analogous 
to  that  of  aqueous  trimethyllead  in  many  respects.  It 
was  determined  in  those  studies  that  an  important  species 
in  the  acid-base  chemistry  of  methylmercury  is  the  dimer 
(CH3Hg)2OH+.  The  asymmetrical  nature  of  the  equilibria 
observed  in  the  methylmercury  system  is  for  the  most  part 
due  to  this  species.  Hence,  a  model  of  the  trimethyllead 
system  including  the  analogous  dimeric  species  ( (CH^) ^Pb) 20h 
was  tested.  This  model,  as  described  by  Equations  9 
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through  12,  gave  the  best  fit  to  the  observed  chemical 
shift  data  of  all  the  models  examined. 

(CH3) jPb+  +  (CH3)3PbOHi  ((CH3)3Pb)  0H+  (11) 

[ ( (CH,) ,Pb) ,0H+] 

K2  =  - | - - -  (12) 

[(CH3)3Pb  ] [(CH3)3PbOH] 

and  are  the  formation  constants  for  the  mononuclear 
and  dinuclear  hydroxide  complexes,  respectively.  No 
evidence  was  found  for  the  existence  of  a  fourth  tri- 
methyllead  species  analogous  to  the  species  (CH^Hg) 
described  for  the  methylmercury  system  (101) . 

The  dinuclear  species  ( (CH^ ) ^Pb) 20H+  was  therefore 
considered  to  form  only  by  the  association  of  (CH^^Pb4^  and 
(CH^J^PbOH  as  represented  by  Equation  11.  Each  trimethyl- 
lead  moiety  in  the  dimeric  species  is  thought  to  share 
a  single  hydroxide  ion  with  the  other;  hence,  the  methyl 
protons  of  the  dimer  would  be  shielded  by  an  amount 
intermediate  to  that  of  the  uncomplexed  (CH^J^Pb  and  the 
singly-complexed  (CH^J-^PbOH.  As  a  result,  the  dinuclear 
species  would  be  expected  to  have  a  chemical  shift  inter¬ 
mediate  to  those  of  the  two  mononuclear  species.  This, 
together  with  the  asymmetric  equilibrium  of  Equation  11, 
fitted  well  with  the  observed  concentration  dependence 
of  the  experimental  data.  The  validity  of  this  model 
was  confirmed  and  and  I<2  were  evaluated  as  well?  by 
comparing  the  observed  chemical  shift  vs  pH  curves  to 
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curves  calculated  from  this  model  by  computer. 

Other  models,  in  which  either  the  species  ((CH3)3Pb)20 
or  ( (CH3) 3PbOH) 2  were  included,  gave  poor  fits.  The 
oxonium  species  ( (CH3 ) 3Pb) 30+  was  not  considered  due  to 
the  low  value  (0.7  M  estimated  for  the  formation 
constant  of  (CH3Hg) 30+  in  the  analogous  methylmercury 
system  (101) ;  the  measurable  formation  constants  of  tri- 
methyllead  complexes  were  found  to  be  consistently  several 
orders  of  magnitude  lower  than  those  of  the  corresponding 
methylmercury  complexes. 

(b)  Determination  of  the  Formation  Constants  of 
Trimethyllead  Hydroxide  Complexes.  The  formation  constants 
and  K as  defined  by  Equations  10  and  12,  and  the 
chemical  shift  of  the  dinuclear  species  ( (CH3 ) 3?b) 30H+ * 

6  ,  were  determined  by  computer  curve-f itting  of  the 

observed  chemical  shift  data  to  the  model  curve  calculated 
as  described  below. 

The  exchange  of  the  trimethyllead  moiety  between  its 
various  aquated  forms  is  fast  on  the  nmr  time  scale; 
hence,  the  observed  chemical  shift  of  its  central  resonance 
(or  its  coupling  constant)  is  an  average  of  the  chemical 
shifts  (or  coupling  constants)  of  the  various  trimethyllead 
species  in  solution,  weighted  according  to  their  individual 
populations.  This  relationship  is  described  by  Equation  1. 
For  the  model  of  the  aqueous  trimethyllead  system  described 
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by  Equations  9  through  12,  Equation  1  becomes 


6obs  ptml6tml  +  ptmloh6tmloh  +  ptmld6tmld  (13 

where  TML,  TMLOH  and  TMLD  refer  to  the  species  (CH^)^Pb+, 
(CH^)^Pt>OH  and  (  (CH^ )  ^Pb)  20H+  /  respectively.  P  refers  to 
the  fraction  of  the  total  trimethyllead  present  in  each 
of  the  forms  and  6  is  the  chemical  shift  of  that  form. 

The  fractions  are  defined  as 


PTML  =  [ (CH3) 3Pb+]/TMLT  (14) 

PTMLOH  =  [(CH3)3PbOH]/™LT  (15) 

PTMLD  =  2[((CH3)3Pb)2OH+]/TMLT  (16) 

where  TMLT  represents  the  total  concentration  of  the  tri¬ 
methyllead  moiety: 

TMLT  =  [ (CH^) 3?b+]  +  [  (CH3) 3PbOH]  +  2 [ ( (CH3 ) 3Pb) 20H+] 

(17) 

The  fractions  PTML0H  and  PTMLD  can  be  expressed  solely  in 
terms  of  K1,  K2 ,  aQH-  and  The  rearrangement  of 

Equation  10  to  Equation  18: 

[(CH3)3PbOH]  =  (CH3)  3Pb+]  (18) 

followed  by  the  division  of  both  sides  of  Equation  18  by 
TMLT  yields  Equation  19. 

P  =  k  a  P  (19) 

1 TMLOH  1  OH"  TML 
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Similarly,  PTMLD  can  be  expressed  as 

PTMLD  =  2KiK2aOH-((CH3)3Pb+)2/™LT  (20) 

,  n 

Equation  14  allows  substitution  of  (TMLT • P  )  for 

TML 

+  2 

((CH^J^Pb  )  to  obtain  Equation  21. 

PTMLD  =  2KlK2aOH“  (™LT)PTML  (21) 

Substitution  of  Equations  19  and  21  into  Equation  13 
gives  Equation  22. 


6OBS  P  TML 6  TML  +  KiaOH“PTML6TMLOH 

+  2KlK2aOH-(™LT)P™L5TMLD  <22> 

PTML  Can  exPressec^  in  terms  of  ,  K 2  and  TMLT  by 
manipulation  of  the  mass  balance  equation  for  the  trimethyl- 
lead  moiety  (Equation  17) .  Dividing  both  sides  of  Equation 
17  by  TMLT  leads  to  Equation  23. 


1 


P  +  P  +  P 

TML  TMLOH  TMLD 


(23) 


Substitution  of  Equations  19  and  21  into  Equation  23, 
followed  by  a  minor  rearrangement,  yields  Equation  24. 

2K1K2aOH-(TMLT)P2ML  +  (1  +  K^oh-^TML-  1  =  0  (24) 

Solution  of  this  equation  for  PTML  by  the  quadratic  formula 
gives  Equation  25. 
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TML 


~(1  +  KlaQH~)  +  ((1  +  KlaQH-)2  +  SK^a^.tTMLT))172 

4KlK2aoH-(TMLT) 


(25) 

The  substitution  of  this  equation  into  Equation  22 

allows  the  evaluation  of  K,  ,  and  from  the  chemical 

1  2  TMLD 

shift  data,  over  the  pH  range  0.5  to  13.0  and  over  the 

-3  -1 

trimethyllead  concentration  range  of  5  x  10  M  to  2  x  10 
M,  using  the  non-linear  least  squares  curve-fitting  program 
KINET  described  by  Dye  and  Nicely  (94) .  The  ionic  strength 
was  controlled  at  0.3  M  with  sodium  perchlorate.  The 
chemical  shifts  of  (CH^J^Pb*  and  (CH^^PbOH  were  taken  to 
be  the  limiting  shifts  of  the  trimethyllead  central  reson¬ 
ance  in  sufficiently  acidic  and  alkaline  solutions, 
respectively. 

The  curve-fitting  procedure  required  Equations  22  and 
25.  These  were  supplied  in  a  subroutine  to  the  curve¬ 
fitting  program  KINET.  Using  reasonable  initial  estimates 
for  K±t  K2  and  ^TMLD '  KINET  calculated  PTML  with  Equation 
25  and  then  predicted  a  value,  ^or  resu^-tant 

chemical  shift  with  Equation  22.  This  procedure  was 
repeated  at  each  data  point  to  predict  a  chemical  shift 


vs.  pH  curve  for  those  particular  values  of  K-^,  K2  and 

6  .  The  observed  chemical  shift  data  set  (i.e. ,  the 

TMLD 

observed  curve)  was  compared  to  the  calculated  data  set 
(i.e.,  the  calculated  curve)  and  the  values  of  K^,  K2  and 
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6TMLD  Were  ad justed  iteratively  by  KINET  to  minimize  the 
sum  of  the  squares  of  the  residuals  (SSR)  over  the  entire 
data  set,  where 


SSR  =  2  (6 6 .  ) 

•  l , CALC  i,OBS; 


(26) 


This  procedure  ultimately  gave  the  values  of  =  7.34 
(±  0.1)  x  104  M_1  and  K2  =  3.12  (±  0.1)  x  101  M_1 
for  the  formation  constants  of  trimethyllead  hydroxide 
and  trimethyllead  dimer,  where  the  uncertainties  are  the 
linear  estimates  of  the  standard  deviation.  The  chemical 
shifts  of  the  species  are  1.535  ppm  ((CH3)3Pb+),  1.253  ppm 


( (CH3 ) sPbOH)  and  1.327  ppm  (  (  (CH3) 3Pb) 2OH+)  ,  quoted  from 
DSS  (i.e.,  the  methyl  resonance  of  sodium  2 , 2-dimethyl- 
2-silapentane-5-sulphonic  acid) . 

Due  to  the  asymmetric  nature  of  Equation  12  (that  is, 
a  dimer  is  formed) ,  the  fractional  population  of  trimethyl¬ 
lead  dimer  in  solution  is  a  function  of  the  total  tri¬ 
methyllead  concentration.  The  species  distribution  as  a 
function  of  pH,  calculated  using  the  above  values  for  K-^ 
and  K2  and  a  value  of  0.100  M  for  TMLT,  is  shown  in 
Figure  3.  The  species  distribution  depicted  by  Figure  4 
was  calculated  in  a  similar  manner,  except  that  TMLT  (the 
total  concentration  of  the  trimethyllead  moiety)  was 
5.00  x  10-3  M.  At  the  latter  concentration , the  maximal 
fractional  concentration  of  the  trimethyllead  moiety  present 
as  (  (CH3) 3Pb) 2OH+  is  0.07,  compared  to  a  value  of  0.46 
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pH 


Figure  3.  pH  dependence  of  the  trimethyllead  species 
distribution  in  a  0.100  M  aqueous  solution,  calculated 
using  and  as  in  the  text. 


PROPORTION 


53 


pH 


Figure  4.  pH  dependence  of  the  trimethyllead  species  distri- 

_3 

bution  in  a  5.00  x  10  M  aqueous  solution,  calculated 


using  and  as  in  the  text. 
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at  the  higher  concentration.  Since  K1  and  K2  were  cal¬ 
culated  from  data  acquired  in  solutions  of  0.3  M  ionic 
strength,  these  curves  must  be  considered  accurate  only 
under  similar  conditions. 

A  very  good  fit  of  the  calculated  data  to  the  observed 

data  was  obtained  with  the  model  described  above.  Thus, 

the  previously  mentioned  analogy  of  the  aqueous  methyl- 

mercury  equilibria  (85,101,102)  to  those  of  trimethyl- 

lead  appears  to  be  consistent  and  apt  in  most  respects. 

As  yet,  however,  few  reports  exist  in  the  literature  for 

comparison  with  the  results  detailed  in  this  chapter. 

Table  1  listed  only  two  other  determinations  of  (in 

both  cases  reported  as  the  K  of  aquated  (CH-. ) -.Pb+)  , 

a  j  j 

both  of  which  can  be  considered  somewhat  unreliable.  The 

4 

first  value  . listed,  3.5  x  10  ,  was  determined  under  condi¬ 
tions  of  unknown  temperature  and  ionic  strength  (32) . 

5 

The  second  value,  5.0  x  10  was  determined  at  25°C  and 
in  0.1  M  KNO^  (72),  but  is  in  doubt  as  well.  The  source 
of  this  reference.  Chemical  Abstracts,  quotes  what  appears 
to  be  the  Ka  of  aquated  trimethyllead  as  8.70.  This  is 
presumably  a  misprint  or  a  misinterpretation  of  the  original 
paper;  it  is  assumed  here  that  this  value  is  actually 

5 

log  K  .  The  value  of  5.0  x  10  for  K±  was  calculated  on 
this  premise  and  is  therefore  questionable.  Nonetheless, 
the  value  determined  for  in  this  chapter  is  inter¬ 
mediate  to  the  two  previously  reported  values.  Other 


than  that  paper  in  which  the  work  described  in  this  paper 
is  outlined  (62)  ,  no  report  of  is  known  to  exist  in 


the  literature. 


CHAPTER  IV 


THE  TRIMETHYLLEAD  COMPLEXES  OF  SELECTED  THIOLS 

AND  AMINO  ACIDS 


A.  Introduction 

This  chapter  presents  the  results  of  a  series  of  experi¬ 
ments  designed  to  investigate  and  quantify  the  binding  of 
trimethyllead  by  certain  organic  ligands.  The  complexation 
of  trimethyllead  by  glycine,  cysteine  and  glutathione 
was  studied  to  partially  elucidate  the  behaviour  of 
trimethyllead  in  physiological  systems,  where  these 
ligands  are  commonly  present.  To  further  increase  our 
understanding  of  sulfhydryl-bound  trimethyllead,  the 
complexes  of  2-mercaptoethanol ,  N-acetylpenicillamine  and 
penicillamine  were  also  considered. 

In  each  section,  the  determinations  of  the  appropriate 
acid  dissociation  constants  are  followed  by  calculations 
of  the  complex  formation  constants  from  chemical  shift  vs. 
pH  measurements.  The  discussion  examines  correlations 
and  trends  in  the  calculated  formation  constants  and 
chemical  shifts,  and  also  presents  species  distribution 
curves  of  the  trimethyllead-complexed  systems. 
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B .  Results 

1 .  2-Mercaptoethanol  Complexes 

(a)  The  Acid  Dissociation  Constant  of  2-Mercapto¬ 
ethanol  .  As  pH  increases,  fully  protonated  2-mercapto- 
ethanol  will  lose  a  proton  as  shown  by  Equation  27. 

HOCH2CH2SH  t  HOCH2CH2S"  +  H+  (27) 

The  acid  dissociation  constant  which  quantifies  this  equili¬ 
brium  is  defined  by  Equation  28. 

a  +  [HOCH-CH  s“] 

K  =  — - - - - -  (28) 

a  [hoch2ch2sh] 

This  is  a  mixed  acid  dissociation  constant:  both  concentrations 

and  activities  are  included  in  the  definition.  All  K  '  s 

a 

used  in  this  thesis  for  calculation  are  of  this  form. 

It  should  also  be  noted  that  at  very  high  pH,  the  species 
HOCH2CH2S~  will  deprotonate  at  the  alcohol  group.  However, 
this  doubly  deprotonated  species  was  found  to  have  no 
significance  in  this  study  because  such  extreme  conditions 
were  never  employed.  For  the  purposes  of  this  work,  the 
acid-base  chemistry  of  2-mercaptoethanol  is  completely 
described  by  Equations  27  and  28. 

The  acid  dissociation  constant,  K^,  was  determined 
by  pH  titration  as  described  in  Chapter  II.  Temperature 
was  maintained  at  25°C  and  ionic  strength  was  controlled 
at  0.3  M  with  NaClC>4.  The  value  of  the  pKa  determined 


from  four  replicates  (corrected  for  hydroxide  ion  bias  as 
recommended  by  Albert  and  Serjeant  (97))  was  9.57  ±  0.02, 
where  the  error  factor  indicates  the  spread  of  the  values 
obtained.  This  was  in  good  agreement  with  the  pK  values 
of  9.5  obtained  by  Calvin  at  25°C  in  0.15  M  NaCl  (103) 
and  9.52  obtained  by  Schwarzenbach  and  Schellenberg  at 
20°C  in  0.1  M  KN03  (102). 

A  value  of  the  pK^  was  also  determined  using  the  nmr 
chemical  shift  vs  pH  method  outlined  in  Chapter  II  and 
detailed  below.  The  observed  chemical  shift  of  the  two 
methylene  groups  of  HOCH^CH^SH  are  the  average  of  the 
corresponding  shifts  of  the  protonated  and  deprotonated 
species,  weighted  according  to  their  populations.  (In  the 
interests  of  brevity,  HOCH2CH2SH  and  HOCH2CH2S  hereafter 
will  be  represented  by  RSH  and  RS  ,  respectively.) 


6OBS 


P  6 
RSH  RSH 


+  P 


(29) 


Both  P^.,,  and  P„„_  will  be  shown  to  be  functions  of  K  ; 
RSH  Rb  a 

hence,  solution  of  Equation  29  enables  the  determination 
of  the  pK^  of  2-mercaptoethanol . 

The  values  of  6RSR  and  6RS_  were  taken  to  be  the 
limiting  chemical  shifts  of  the  appropriate  resonances 
at  low  and  high  pH,  respectively.  The  fractional  popu¬ 
lations  are  defined  as 


p 

RSH 


TOT 


[ RSH ] /  [ R] 


(30) 


* 


[RS  ]/[R] 


TOT 


(31) 


PRS“ 


where  [R]^^  is  the  total  concentration  of  2-mercapto- 
ethanol  in  all  forms.  The  summation  of  the  fractional 


populations  is  unity: 


1 


P  +  P 
RSH  RS" 


(32) 


Appropriate  substitution  of  Equations  30  and  31  into 
Equation  28  results  in  an  expression  of  in  terms  of  the 
fractional  populations  defined  above. 


(aH+*PRS-)//  PRSH 


(33) 


Appropriate  substitution  and  rearrangement  of  Equations  32 

and  33  gives  expressions  of  the  fractional  populations  in 

terms  of  hydrogen  ion  activity  and  K  alone. 

a 


RSH  aH+^(Ka  +  aH+* 

(34) 

RS-  *  V(Ka  +  aH+> 

(35) 

Substitution  of  these  equations  into  Equation  29  results 

in  an  equation  for  the  observed  chemical  shift  in  terms 

of  known  experimental  parameters  and  only  one  unknown, 

K  . 
a 


OBS 


=  (aH+ 


6RSH  +  Ka 


{RS-)/(Ka  +  aH+) 


(36) 


With  this  equation,  pK_^  determined  from  chemical  shift 

vs.  pH  titration  data. 
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The  pmr  spectrum  of  2-mercaptoethanol  consists  of 
two  multiplets  of  approximately  equal  intensity.  The 
methylene  protons  nearest  the  alcohol  group  (-0CH2-)  give 
rise  to  the  multiplet  further  upfield,  while  those  nearest 
the  sulfhydryl  (-CH2S-)  cause  the  downfield  multiplet. 
Figure  5  shows  the  pH  dependence  of  these  two  resonances. 

K  was  evaluated  with  Equation  36  using  the  non-linear 
least-squares  curve-fitting  program  KINET  (94)  described 
earlier  in  this  thesis.  Each  curve  in  Figure  5  was  con¬ 
sidered  separately,  over  the  pH  range  5.9  to  13.0.  Total 
2-mercaptoethanol  concentration  was  0.139  M,  temperature 
was  25°C  and  ionic  strength  was  maintained  at  0.3  M  with 
NaClO^.  Using  the  limiting  chemical  shift  values  of 
3.706  ppm  vs.  DSS  for  6RSR  and  3.537  ppm  vs.  DSS  for  6RS_, 
the  -CH2S-  data  gave  a  value  for  Ka  of  2.47  (±  0.04)  x  10  ^ 
M.  Similarly,  using  2.683  ppm  vs.  DSS  for  6RRR  and  2.566 
ppm  vs.  DSS  for  SRS_,  the  -OCH2"  data  resulted  in  a  Ka 
of  2.31  (±  0.04)  x  10-10  M.  The  averaged  result  gave  the 
final  pKa  of  9.62  ±  0.02.  The  conditions  and  methods 
employed  in  the  determination  closely  match  those  used 
in  the  determination  of  the  formation  constants  of  the 
2-mercaptoethanol- trimethyl lead  complexes  described  later 
in  this  section;  therefore,  the  value  used  for  the  pKa  of 
2-mercaptoethanol  in  this  study  was  9.62  ±  0.02  rather  than 
the  value  of  9.57  obtained  by  pH  titration. 


to 


CHEMICAL  SHIFT,  ppm  vs  DSS 
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PH 


Figure  5.  pH  dependence  of  the  chemical  shift  of  the 
-OCH2-  and  the  -CH2S-  multiplet  resonances  for  a 
0.139  M  solution  of  2-mercaptoethanol  at  25°C  and  0.3 
M  ionic  strength  (maintained  with  NaClO^) . 
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(b)  The  2-Mercaptoethanol  Complexes  of  Trimethyllead . 
The  binding  of  trimethyllead  by  2-mercaptoethanol  was 
investigated  by  monitoring  the  chemical  shift  of  the  methyl 
protons  of  trimethyllead  as  a  function  of  pH  and  also  as 
a  function  of  the  trimethyllead  to  2-mercaptoethanol 
concentration  ratio. 

(i)  Qualitative  Observations 

The  chemical  shift  of  the  methyl  protons  of  trimethyl¬ 
lead  is  presented  in  Figure  6  as  a  function  of  pH  for  three 
solutions  containing  trimethyllead  and  2-mercaptoethanol 
and  for  one  solution  of  trimethyllead  alone.  The  curves 
for  the  solutions  containing  2-mercaptoethanol  approach 
the  curve  for  the  solution  containing  only  trimethyllead 
at  the  extremes  of  pH.  These  results  indicate  that,  at 
low  pH,  protons  compete  with  the  trimethyllead  for  the 
ligand  sulfhydryl  group,  while  at  high  pH,  hydroxide  ions 
compete  effectively  with  the  ligand  for  the  available 
trimethyllead  cations.  The  attainment  of  maximum  complexa- 
tion  is  indicated  by  the  plateaux  in  the  intermediate  pH 
range,  where  both  competing  ions  are  at  low  concentrations. 

A  mole  ratio  experiment  was  performed  to  determine 
the  coordination  numbers  of  these  complexes.  The  results 
of  this  experiment,  in  which  the  chemical  shift  of  tri¬ 
methyllead  was  measured  as  a  function  of  the  ligand  to 
metal  ratio,  are  shown  in  Figure  7.  The  observed  chemical 
shift  varies  with  increasing  proportion  of  ligand  until  a 
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Figure  6.  pH  dependence  of  the  chemical  shift  of  trimethyl- 
lead  in  a  sulfhydryl  —  free  0.0100  M  aqueous  solution  (A) 
and  in  solutions  of  0.0100  M  trimethyllead  and  0.0200  M 
2-mercaptoethanol  (•),  0.0100  M  trimethyllead  and  0.00970 
2-mercaptoethanol  (♦),  and  0.0200  M  trimethyllead  and 
0.0100  M  2-mercaptoethanol  (■).  Ionic  strength  was  con¬ 
trolled  at  0.3  M  with  NaClC>4  and  temperature  was  25°C. 
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Figure  7.  The  dependence  of  the  chemical  shift  of  tri- 
methyllead  upon  the  ratio  of  2-mercaptoethanol  to  tri- 
methyllead  in  solution  at  pH  8.0.  The  concentration  of 
ligand  was  0.00887  M  except  at  the  ratio  of  zero  to  one. 
The  ionic  strength  was  maintained  at  0.3  M  with  NaClO^  and 
the  temperature  was  25°C. 
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one-to-one  ratio  is  reached,  after  which  the  chemical  shift 
remains  constant  even  with  a  several— fold  excess  of  ligand. 
This  indicates  that  the  major  complex  at  all  measured 
ratios  is  the  one— to— one  ligand-to-metal  complex, 

(CH3 )  sPbSCF^Cf^OH .  No  evidence  was  found  for  any  complex 
of  greater  ligand  to  metal  ratio.  However,  the  deviation 
from  linearity  in  the  low  ligand  to  metal  portion  of  the 
curve  was  attributed  to  the  presence  of  a  weak  1  to  2 
ligand  to  metal  complex,  most  likely  ( (CH3) 3Pb) 2SR+.  This 
was  confirmed  by  studies  described  later  in  this  section. 

(ii)  derivation  of  a  Mathematical  Model  for  the 
2-Mer cap  toe thanol -Trim ethy Head  System 

Separate  pmr  signals  are  not  observed  for  the  different 
trimethyllead  species,  indicating  that  exchange  of  the 
trimethyllead  moiety  among  its  various  forms  is  fast  on 
the  nmr  time  scale  (62)  .  The  chemical  shift  of  the  observed 
resonance  is  an  average  of  the  chemical  shifts  of  the 
various  trimethyllead  species  present  in  solution,  weighted 
according  to  their  populations.  Formation  constants  and 
chemical  shifts  for  the  metal-ligand  complexes  were 
calculated  by  fitting  proposed  models  of  the  complexation 
equilibria  to  the  observed  chemical  shift  data,  using 
the  curve— f itting  computer  program  KINET  (94)  described 
earlier  in  this  thesis.  The  best  fit  to  the  observed 
data  was  obtained  for  a  model  which  included  one  to  one 
and  two  to  one  trimethyllead  to  2-mercaptoethanol  complexes. 
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The  aim  in  this  section  of  the  thesis  is  to  derive  a 
mathematical  model  for  the  observed  trimethyllead  chemical 
shift  behaviour;  this  model  will  be  used  in  the  next  section 
to  calculate  the  unknown  parameters  for  the  metal-ligand 
complexes . 

The  derivation  of  the  model  of  any  chemical  system 
requires  the  initial  definition  of  equations  which  quantify 
the  component  equilibria  of  the  system.  The  equilibria 
of  trimethyllead  in  aqueous  solutions  are  described  by 
Equations  9  through  12  in  Chapter  III;  those  of  2-mercapto- 
ethanol  are  given  by  Equations  27  and  28  in  this  chapter. 
Equations  37  through  40  describe  the  equilibria  of  the 
trimethyllead-2-mercaptoethanol  complexes  which  resulted 
in  the  model  most  closely  fitting  the  observed  chemical 
shift  data. ■ 

TML+  +  RS~  t  TMLSR  (37) 


K 


fS 


[TMLSR] 
[TML+] [RS~] 


(38) 


TML+  +  TMLSR  t  (TML) 2SR+ 

[ (TML) 2sr+] 

—  _ — - - 

fD  [TML+] [TMLSR] 


(39) 

(40) 


TMLSR  and  (TML) 2SR+  refer  to  (CH3) 3PbSCH2CH2OH  and 
( (CH3) 3Pb) 2SCH2CH2OH+,  respectively,  and  Kfg  and  KfD 


are 


< 
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the  formation  constants  of  these  complexes.  The  exchange- 
averaged  chemical  shift  of  this  system  is  given  by  Equation 
41;  this  equation  is  the  fundamental  chemical  shift 
relation  of  this  model. 


^OBS  P TML ^ TML  +  PTMLOH^TMLOH  +  PTMLD^TMLD 

+  PTMLSR6TMLSR  +  P (TML) 2SR+6 (TML) 2SR+  (41) 

The  first  three  terms  are  defined  as  for  Equations  13 
through  16  in  Chapter  III.  The  fractions  of  the  remaining 
species  are  defined  as  follows: 


TMLSR  =  [™LSR]/TMLT 

(42) 

(TML)2SR+  =  2[(TML)2SR+1/TMLT 

(43) 

where  TMLT  is  the  total  concentration  of  the  trimethyllead 
moiety  in  solution.  Appropriate  substitution  into  Equation 
41  results  in  Equation  44. 


6OBS  PTML6TML  +  KiaOH-PTML6TMLOH 

+  2K1K2a0H_(™LT^  PTML6TMLD  +  KfS[RS  ]PTML  6TMLSR 
+  2KfDKfS[RS  1 (TMET) pTML  6 (TML) 2SR+  (44) 

where 

P  =  [TML+]/TMLT  (45) 

TML  1 

All  that  remains  is  to  solve  for  [TML+]  and  [RS  ] 


* 
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in  terms  of  known  parameters.  The  mass  balance  equations 
for  total  trimethyl lead  and  total  2-mercaptoethanol  provide 
the  necessary  relations: 

TMLT  =  [TML+]  +  [TMLOH]  +  2 [TMLD]  +  [TMLSR] 

+  2 [ (TML) 2SR+]  (46) 

RSHT  =  [RS~]  +  [RSH]  +  [TMLSR]  +  [  (TML)  2SR+]  (47) 


where  RSHT  is  the  total  concentration  of  the  2-mercapto- 
ethanol  moiety. 

By  appropriate  substitution  of  Equations  10,  12,  28, 
38  and  40,  the  mass  balance  equations  can  be  written  as 
two  equations  in  two  unknowns,  [TML+]  and  [RS  ] .  These 
can  be  manipulated  to  give  a  polynomial  in  [TML+]  only, 
the  solution  of  which  allows  the  determination  of  both 
[TML+]  and  [RS-]  if  TMLT,  RSHT  and  pH  are  known.  This 


derivation  is  outlined  below. 


The  mass  balance  equation  for  the  ligand,  expressed  as 


+ 


a  function  of  [TML  ]  and  [RS  ]  only,  is  equation  48. 


RSHT  =  [RS  ]  + 


aH+ [Rs  ] 


K 


+  Kf s [TML  ] [RS  ] 


+  KfDKfs(TML+]2[RS  ] 


(48) 


Rearrangement  of  this  gives  Equation  49: 

[RS~]  =  RSHT/  ( 1  +  aH+//Ka  +  KfS[™L+] 
+  KfDKfs[TML+]2) 


(49) 


4 


. 


Similar  substitution  into  the  mass  balance  equation  for 
trimethyllead  (Equation  46)  results  in  Equation  50: 


TMLT  =  [TML+]  +  K1a0H_[TML+]  +  2K1K2a  [TML+] 2 

+  Kfs(TML+] [RS-]  +  2KfDKfs[TML+]2[RS_]  (50) 

Substitution  of  the  right  hand  side  of  Equation  49  for 
[RS  ]  in  Equation  50  gives  an  expression  for  the  mass 
balance  of  trimethyllead  in  terms  of  [TML+]  alone. 

TMLT  =  [TML+]  +  K1aQH_ [TML+]  +  2K1K2aQH_ [TML+] 2 

+  KfstTML+] {RSHT/ ( 1  +  aH+/Ka  +  Kf g [TML+] 

+  KfDKfs[TML+]2)}  +  2KfDKfs[TML+]2{RSHT/ 

(1  +  aH+/Ka  +  Kfs[TML+)  +  KfDKfs(TML+]2)}  (51) 

A  rearrangement  of  this  expression  gives  a  polynomial  in 
decreasing  powers  of  [TML+] : 

{2K1K2aOH_KfDKfs} [TML  ]  +  ( 2KxK2aOR-Kf S 

+  (1  +  Kxa0H~*  KfDKfS^™L  '  +  '2KlK2aOH"'(1  +  aH+//Ra) 

+  (1  +  KiaoH_)Kfs  +  2KfDKfs(RSHT)  -  Kf[)Kf £  (TMLT)  }  [TML+] 

+  {(1  +  K^.1  (1  +  aH+/Ka)  +  Kfs  (RSHT)  -  Kf£  (TMLT)  } 
[TML+]  -  {(1  +  aR+/Ka) -(TMLT)  )  =  0  (52) 


. 
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The  coefficients  of  this  polynomial  are  functions  of  K^, 

K9,  K  (known  constants),  RSHT ,  TMLT ,  a„+  and  a^„_  (known 
experimental  parameters)  and  the  formation  constants  Kfg 
and  the  values  of  which  were  initially  estimated 

but  later  optimized  by  computer  calculations.  Given  this 
information.  Equation  52  can  be  solved  for  [TML+] ;  hence, 
the  concentration  of  all  the  above-mentioned  metal  and 
ligand  species  can  be  determined.  Finally,  the  substitution 
of  the  calculated  values  of  [TML+]  from  Equation  52  and 
[RS  ]  from  Equation  49  allows  the  calculation  of  the 
observed  chemical  shift  for  any  reasonable  set  of  solution 
conditions,  including  those  given  in  Figure  6. 

(Hi)  Computer  Optimization  of  Formation 
Constants  and  Chemical  Shifts  of  Trimethyllead  Complexes 

Due  to  the  complexity  of  the  equations  derived  above, 

the  determinations  of  Kfg,  Kf[),  «TMLSR  and  5  (TML)  2sr+ 
were  accomplished  with  the  digital  computer  program  KINET 
(94) ,  described  earlier  in  this  thesis.  A  flow  chart 
summarizing  this  optimization  procedure  is  shown  in 
Figure  8.  Initially,  the  program  was  supplied  the  known 

values  of  K2 ,  Kg,  «TML/  5TmlOH  and  5TMLD'  ExPeriments 

carried  out  at  various  metal  to  ligand  ratios  (as  plotted 

in  Figure  6)  provided  the  required  chemical  shift  vs.  pH, 
TMLT  and  RSHT  data.  In  addition,  initial  estimates  were 
introduced  for  the  formation  constants  and  limiting 
chemical  shifts  of  the  trimethyllead  complexes. 
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Figure  8.  Flow  chart  of  the  computer  procedure  employed 
to  optimize  the  formation  constants  and 
chemical  shifts  of  the  2-mercaptoethanol 


complexes  of  trimethyllead 
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The  mathematical  model  (Equation  44)  ,  along  with  all 
preliminary  calculations,  was  supplied  to  KINET  in  a  sub¬ 
routine  (hereafter  called  EQN) .  At  each  data  point,  the 
subroutine  EQN  calculated  the  coefficients  of  Equation 
52  and  supplied  them  to  the  IBM  SSP  Library  Subroutine 
DRTMI ,  which  solved  the  polynomial  for  a  value  of  [TML+] 
between  zero  and  TMLT.  This  enabled  the  successive  solution 
of  Equations  45,  49  and  44,  thereby  providing  a  value  for 
the  averaged  chemical  shift,  <$CALC,  °f  the  trimethyllead 
protons  for  one  data  point  and  a  particular  set  of  esti¬ 
mated  parameters.  The  residual,  which  is  the  difference 
between  this  calculated  shift  and  the  experimentally 
observed  shift,  was  determined  for  these  conditions. 

This  process  was  repeated  until  the  residual  at  every  data 
point  had  been  calculated.  KINET  then  computed  the  sum 
of  the  squares  of  the  residuals  (SSR)  as  defined  by 
Equation  26  in  Chapter  III.  After  each  such  calculation, 
KINET  altered  the  values  of  the  estimated  parameters 
such  that  SSR  would  be  less  in  the  following  calculation. 

The  best  values  for  Kfg,  KfD,  «TMLSR  and  6 (tml) 2SR+  W6re 
attained  when  the  minimum  in  the  SSR  was  reached.  At 
this  time,  KINET  also  provided  error  estimates  for  these 
parameters  based  on  the  quality  of  the  fit  between  the 
observed  and  calculated  chemical  shift  curves. 
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(iv)  Computed  Results 

The  value  of  Kf  was  computed  to  be  9.05  (±  0.11)  x  105, 
using  all  the  data  in  Figure  6  over  the  pH  range  1.5  to 
12.5  simultaneously.  Similarly,  the  best  values  of  Kf  , 

6TMLSR  and  6(TML)2SR+  Were  calculated  to  be  14.0  ±  8.8, 
1.3245  i  0.0002  and  1.454  ±  0.007  respectively,  where  the 

chemical  shifts  are  in  units  of  ppm  vs.  DSS.  The  uncer¬ 
tainties  are  linear  estimates  of  the  standard  deviations 
as  determined  by  KINET. 

No  other  determinations  of  these  quantities  are 
available  in  the  literature  for  comparison.  However,  the 
good  quality  of  the  fit  of  this  model  calculation  to  the 
observed  data  tends  to  support  the  accuracy  of  the  optimized 
parameters,  as  the  calculated  error  estimates  are  a  fair 
measure  of  the  precision  and  accuracy  of  these  results. 

All  are  of  reasonably  low  magnitude  except  for  that  of 
Kf d /  this  large  relative  error  is  due  in  part  to  the  low 
magnitude  of  KfD-  The  concentration  of  (TML)2SR+  is 
always  a  small  fraction  of  the  total  concentration  and  con¬ 
sequently  cannot  be  determined  with  a  high  degree  of 

precision.  The  fraction  of  (TML) 2SR+  is  reduced  further 
by  the  unfavorable  ratios  of  metal  to  ligand  employed  for 


some  of  the  experimental  solutions. 

No  evidence  was  found  to  indicate  the  presence  of 
measurable  quantities  of  any  species  in  solution  other  than 

those  described  above. 
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2 .  Glycine  Complexes 


The  binding  of  trimethyllead  by  glycine  was  studied 
as  a  model  for  complexation  by  the  amino  and  carboxylate 
groups  in  the  sulf hydryl-containing  amino  acids  and  pep¬ 
tides  studied  later  in  this  thesis. 

(a)  The  Acid  Dissociation  Constants  of  Glycine. 
Fully  protonated  glycine  is  a  dibasic  acid,  dissociating 
as  follows: 


H2NCH2COOH 

-> 

H^NCH2COO 

+  H+ 

(53) 

H2NCH2COO~ 

-> 

H2NCH2COO" 

+  H+ 

(54) 

The  corresponding  acid  dissociation  constants  are  defined 
below: 


[HL]aR+ 

[h2l+] 

[L  ] aH+ 
[HL] 


(55) 

(56) 


where  H2L+,  HL  and  L  represent  the  fully  protonated, 
zwitterionic  and  doubly  deprotonated  species  respectively. 
The  numerical  values  for  and  K^2  were  determined  by 

monitoring  the  chemical  shift  of  the  glycine  methylene 
protons  as  a  function  of  pH  at  25°C  in  a  0.100  M  glycine 
solution,  controlled  at  0.3  M  ionic  strength  with  NaC104- 
The  resultant  curve  is  plotted  in  Figure  9. 
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pH 


Figure  9.  pH  dependence  of  the  chemical  shift  of  the 
methylene  protons  of  glycine  in  a  0.100  M  solution  at  25°C 
and  0.3  M  ionic  strength  (maintained  with  NaC104) . 
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The  mathematical  model  of  glycine  acid  dissociation 
was  derived  in  a  manner  quite  analogous  to  that  of  2- 
mercaptoethanol ,  except  that  two  acidic  protons  were 
considered  rather  than  just  one.  The  values  of  pKa^  and 
pKa2  are  sufficiently  separated  that  the  first  ionization 
is  essentially  complete  before  the  second  ionization 
begins  to  occur.  This  enables  the  calculation  of  both 
these  constants  as  simple  independent  ionizations  in  the 
same  run  of  the  curve-fitting  computer  program  KINET. 

Using  chemical  shifts  of  3.552  ppm  vs.  DSS  for  zwitterionic 
glycine  and  3.175  ppm  vs.  DSS  for  doubly  deprotonated 
glycine,  values  of  3.44  (±  0.06)  x  10  ^  and  1.90  (f  0.01) 
x  10~10  were  obtained  for  K  2  and  Ka2,  respectively,  and 
a  value  of  3.906  ±  0.001  ppm  vs.  DSS  was  obtained  for  the 
chemical  shift  of  fully  protonated  glycine.  From  these 
values,  pKa-^  oalculated  to  2.46  and  pKa2 
9.72,  which  compare  favorably  with  the  recently-published 

(104)  pKal  of  2.34  and  pKa2  of  9.69,  determined  under 
conditions  of  0.2  to  0.3  M  ionic  strength  and  25 °C  tempera¬ 
ture.  Other  values  of  these  constants  found  in  the 
literature  include  a  pK^  of  2.45  and  pKa2  °f  9.75, 
determined  spectrophotometrically  at  25  C  in  1  M  NaClO^ 

(105) ;  and  a  pKal  of  2.47  and  pKa2  of  9.68,  determined 
using  a  glass  electrode  at  25°C  in  0.2  M  KC1  (106) .  The 
values  determined  by  this  pmr  technique  were  used  in  the 
following  computation  of  the  formation  constants  of 
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the  glycine  carboxyl  complexes  of  trimethyllead. 

(b)  The  Glycine  Complexes  of  Trimethyllead . 

The  formation  constants  and  chemical  shifts  of  the  glycine 
complexes  of  trimethyllead  were  determined  by  measuring 
the  chemical  shift  of  the  methyl  protons  of  trimethyllead 
as  a  function  of  pH.  The  chemical  shift  vs.  pH  curve 
resulting  from  a  solution  having  a  metal  to  ligand  ratio 
of  one  to  ten  is  compared  to  that  of  a  simple  trimethyllead 
solution  in  Figure  10.  Ratios  other  than  one  to  ten  were 
attempted,  using  both  the  trimethyllead  methyl  and  the 
glycine  methylene  proton  chemical  shifts  as  functions  of 
pH.  However,  due  to  the  low  degree  of  complexation  of 
trimethyllead  by  glycine,  only  the  most  extreme  ratio  of 
metal  to  ligand  resulted  in  a  curve  measurably  different 
from  that  for  a  solution  free  of  metal-ligand  complex. 
Therefore,  the  program  KINET  could  calculate  formation 
constants  only  with  the  data  shown  in  Figure  10. 

The  following  reaction  scheme  was  found  to  give  the 
best  fit  to  the  observed  chemical  shift  data: 

Kai  Ka2  _ 

hInch.cooh  i  h+nch9cocT  t  H?NCH2COO 

3  2  •  3  2  r,  „  <57> 

H  Kfl  Ka3  +  +  Kf 2 
HtNCH2C00Pb (CH3) 3  t  H2NCH2COOPb(CH3) 3 

where  Kfl  and  Kf2  are  the  formation  constants  of  the 
complexes  of  trimethyllead  with  monoprotonated  and  fully 
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Figure  10.  pH  dependence  of  the  chemical  shift  of  tri¬ 
methyllead  at  25°C  and  0.3  M  ionic  strength  (maintained 
with  NaC104)  in  a  solution  of  0.0100  M  trimethyllead  only 
(A)  and  in  a  solution  containing  0.0100  M  trimethyllead 
and  0.100  M  glycine  (♦). 
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deprotonated  glycine,  respectively.  The  binding  site  in 
both  cases  is  assumed  to  be  the  carboxyl  group  (as  discussed 
in  part  C  of  this  chapter) .  Using  the  aforementioned 
abbreviations : 


[TML • LH+] 
[TML+] [HL] 
[TML-L] 
[TML+] [L“] 
[TML-L] aR+ 
[TML • LH+] 


(58) 


(59) 


(60) 


The  formation  constants  and  Kf2/  as  well  as  the 

associated  chemical  shifts,  were  determined  by  a  compu¬ 
tation  quite  similar  to  that  described  for  the  2-mercapto- 
ethanol  complexes  of  trimethyllead .  Using  the  appropriate 
constants  detailed  earlier  in  this  thesis,  Equation  61 
was  solved  by  KINET  for  the  optimum  values  of  these 
formation  constants  and  chemical  shifts. 


^OBS  ^TML^TML  +  PTMLOH^TMLOH  +  PTMLD^TMLD 

+  ptml.lh+6tml.lh+  +  ptml.l6tml.l  (61) 

The  values  which  resulted  in  the  best  fit  were  1.26  ±  0.58 
for  Kfl,  1.481  ±  0.020  ppm  vs.  DSS  for  <5TML#LH+,  29*5  ±  3-8 
for  Kf 2  and  1.433  ±  0.016  ppm  vs.  DSS  for  5TML>L. 

No  previously  determined  values  for  these  constants 
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were  available  for  comparison,  but  the  fit  to  the  observed 
data  is  reasonably  good.  The  relative  errors  estimated  here 
by  KINET  are  larger  than  those  found  for  the  2-mercapto- 
ethanol  results,  but  this  is  because  the  competition  by 
protons  and  hydroxide  ions  precludes  the  formation  of  the 
weak  glycine  complexes  through  most  of  the  experimental 
pH  range.  All  the  glycine  results  are  known  with  a  lesser 
degree  of  confidence  due  to  the  relative  scarcity  of 
complex,  which  in  turn  causes  a  higher  error  in  the 
chemical  shift  data  from  which  the  results  must  be  cal¬ 
culated.  Nonetheless,  this  model  still  provides  the  best 
fit  of  the  various  types  investigated. 

The  small  size  of  the  formation  constant  Kf2  is 

of  immediate  practical  use.  This  constant  is  negligible 

5 

compared  to  the  K^g  (9.05  x  10  )  of  the  2-mercaptoethanol 
sulfhydryl  complex  of  trimethyllead .  In  sulfhydryl- 
containing  amino  acids  and  peptides,  the  amino  and  the 
sulfhydryl  groups  ionize  more  or  less  simultaneously; 
assuming  a  similar  large  difference  in  K^g  and  the 

effect  of  the  latter  may  be  neglected  in  all  solutions  of 
low  metal  to  ligand  ratios.  In  addition,  no  evidence 
of  amino  group  complexation  has  been  observed;  such 
complexes  may  therefore  be  neglected  as  well.  These 
assumptions  result  in  a  considerable  simplification  of 
the  computer  calculations  in  the  following  sections, 
without  decreasing  the  accuracy  of  the  determinations  of 
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Kfs  values  and  the  associated  chemical  shifts  for  the 

trimethyllead-amino  acid  complexes.  These  simplifications 

were  employed  as  required  in  the  remainder  of  this  chapter. 

-9 

The  value  of  K  0  was  calculated  to  be  4.4  x  10 

ai 

(pK  =  8.35),  using  the  relation  which  states  that  the 
a  j 

product  of  and  K  ^  equals  the  product  of  K ^  and  K 

This  relation  results  from  the  properties  of  the  cyclic 
system  (118)  and  is  hence  exactly  analogous  to  Equation 
91,  described  and  justified  in  detail  in  section  4(b)  of 
this  chapter. 


3.  N-Acetyl-D ,L-Penicillamine  Complexes 
(a)  The  Acid  Dissociation  Constants  of  N-Acetyl-D , L- 
Penicillamine.  The  structural  formula  of  fully  protonated 
N-acetylpenicillamine  is  given  by  I: 

O  0 

II  II 

H-jCCNHCHCOH  I 

3  I 

C(CH3)2SH 

hereafter  represented  by  t^A  or  by  a  shorthand  version  of 
the  structure  showing  only  the  proton  and  metal  binding 
sites.  As  pH  is  increased,  first  the  carboxyl  and  then 
the  sulf hydryl  protons  are  removed.  The  equilibria  involved 
are  described  by  Equations  62  and  63;  the  corresponding 
acid  dissociation  constants  are  defined  by  Equations  64 


and  65. 
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i— COOH 
SH 


(h2a) 


-COO 
— SH 
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+  H 


[ 


K 


K 


] aH+ 


a2 


[HA  ] 


COO 

—COO 

SH 

— s“ 

(HA_) 

(A= 

[HA-]aH+ 

al 

[H.A] 

+  H 


(62) 


(63) 


(64) 

(65) 


The  values  of  K  ^  and  K^2  were  determined  by  pH  titra¬ 
tion  as  described  in  Chapter  II.  The  ionic  strength  of 

the  solution  was  maintained  at  0.3  M  with  NaClO.  and  tem- 

4 

perature  was  controlled  at  25°C.  The  values  obtained  for 
and  K  2  were  corrected  for  the  interference  of  protons 
and  hydroxide  ions,  respectively,  using  the  methods  of 
calculation  detailed. by  Albert  and  Serjeant  (97).  Four 
replicates  gave  numerical  values  of  3.48  -  0.01  for  pK  ^ 
and  10.28  t  0.05  for  pKa2,  where  the  uncertainty  is  the 
maximum  deviation  from  the  average.  Wilson  and  Martin 
(107)  obtained  the  pK  values  of  this  ligand  by  a  similar 
pH  titration.  In  0.15  M  KC1  at  25°C,  they  determined  a 
pK  _  of  3.3  and  a  pK  0  of  9.9;  these  compare  reasonably 
well  with  the  results  given  above,  considering  the  dif— 
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f erences  in  ionic  strength  and  in  calculations.  There  is 
no  indication  that  these  authors  performed  the  recommended 
(97)  corrections  for  proton  and  hydroxide  ion  concen¬ 
tration  in  their  pK  calculations.  Furthermore,  the  use  of 
only  two  significant  figures  probably  indicates  their 
level  of  confidence  in  the  accuracy  of  the  results.  It 
w.as  confirmed,  however,  that  no  evidence  exists  for  amide 
hydrogen  ionization  below  pH  11  (107);  such  a  possibility 
is  neglected  henceforth. 

The  determination  of  pKa  values  using  nmr  titrations 
of  pH  vs.  chemical  shift  was  not  successful  due  to  the 
low  solubility  of  the  ligand  under  the  required  experimental 
conditions,  combined  with  the  limited  usefulness  of  the 
broad,  low  amplitude  nmr  resonances  of  the  ligand.  This 
lack  of  sharpness  is  most  likely  related  to  the  presence 
of  the  asymmetric  center  near  the  methyl  groups  of  the 
ligand  penicillamine  residue:  the  low  broad  resonances 
are  probably  unresolved  multiplets  (108).  The  combination 
of  poor  spectral  and  solubility  characteristics  was  such 
that  the  accurate  measurement  of  the  chemical  shifts  of 
the  ligand  resonances  as  a  function  of  pH  was  not  possible 
with  the  nmr  instrument  used  in  this  research. 

(b)  The  N-Acetyl-D , L-Penicillamine  Complexes  of 
Tr imethy Head .  The  pH  vs.  chemical  shift  behaviour  of 


the  methyl  protons  of  trimethylload  in  solutions  of  excess 
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N—acety lpenici llamine  was  studied  to  determine  the  type 
and  the  extent  of  metal-ligand  complexation  in  this  system. 
Excess  ligand  was  employed  to  ensure  accurate  results 
for  the  strong  complexes  with  minimal  interference  from 
any  weak  complexes  occurring  in  the  same  pH  range.  Low 
experimental  signal  amplitudes  precluded  the  study  of  the 
pH  vs.  chemical  shift  behaviour  of  the  various  N-acetyl- 
penicillamine  protons. 

The  results  of  experiments  at  metal  to  ligand  ratios 
of  one  to  two  and  one  to  four  are  presented  in  Figure  11. 

The  plateaux  between  pH  8  and  11  indicate  the  formation 
of  a  strong  complex  of  shift  intermediate  to  those  of  the 
uncomplexed  and  the  hydroxyl  complexed  trimethyllead 
cations.  The  difference  in  the  breadth  of  these  plateaux 
is  due  to  the  differences  in  total  ligand  concentration: 
a  higher  concentration  of  ligand  drives  the  complexation 
equilibria  towards  increased  formation  of  metal-ligand 
complexes  over  a  wider  pH  region,  despite  the  competition 
of  protons  and  hydroxyl  ions  mentioned  previously. 

The  model  found  to  generate  the  best  fit  with  the 
observed  chemical  shift  data  was  derived  using  the  following 


system  of  equilibria 
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1  3  5  7  9  11  13 


pH 


Figure  11.  pH  dependence  of  the  chemical  shift  of  tri- 
methyllead  at  25°C  and  0.3  M  ionic  strength  (maintained 
with  NaCl04)  in  solutions  of  0.0100  M  trimethyllead  only 
(A),  0.0100  M  trimethyllead  and  0.0198  M  N-acetyl-D , L- 
penicillamine  (♦),  and  0.0100  M  trimethyllead  and  0.0395 

M  N-acetyl-D / L-penicillamine  (  I)  • 
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where  the  species  TML.AH  and  TML.A  are  the  only  two 
trimethyllead-N-acetylpenicillamine  complexes  present  in 
measurable  amounts.  The  formation  constants  and  K. 

LU  1  u 

are  defined  in  Equations  67  and  68: 


K 


fO 


K 


fS 


[TML.AH] 
[TML+] [HA_] 
[TML.A-] 
[TML+] [A= ] 


(67) 


(68) 


No  higher  order  complexes  containing  more  than  one  metal 
or  ligand  moiety  per  molecule  of  complex  were  detected. 
However,  at  different  concentration  ratios,  this  model  may 
require  modification  to  include  possibilities  such  as 
(TML) 2A. 

The  derivation  and  calculation  of  the  model  was  per¬ 
formed  quite  analogously  to  that  described  for  2-mercapto- 
ethanol  earlier  in  this  chapter.  The  basic  relation 
describing  the  observed  chemical  shift  behaviour  is 
Equation  69: 


^OBS  PTML^TML  +  PTMLOH^TMLOH  +  PTMLD^TMLD 


r* 

+  PTML.AH5TML. AH  +  TML.A-  TML.A' 


(69) 


Using  the  values  of  ^.nd  K^2  calcul3ted  the  previous 

subsection,  the  computer  program  KINET  solved  this  equation 
for  the  best-fitting  values  of  K^q/  and  the  trimethyl- 

lead  pmr  chemical  shifts  of  the  corresponding  complexes, 


87 


i •  e .  ,  omMT  and  6mXKT  These  are  4.0  t  2.1  for  K 


TML.AH  TML .A” * 

1.465  t  0.030  ppm  vs.  DSS  for  6 


fO' 


TML . AH , 


4.07  (±  0.21)  x  10 


for  K^_  and  1.351  ±  0.0005  ppm  vs.  DSS  for  6m..,  ,  _  ,  where  the 
uncertainties  are  linear  estimates  of  the  standard  devia¬ 
tion  provided  by  KINET.  The  larger  relative  errors  cal¬ 
culated  for  K,-_  and  6m..T  are  functions  of  the  scarcity 
of  complex  giving  rise  to  the  original  data.  Nonetheless, 
some  carboxyl-complexed  trimethyllead  must  be  present 
because  omission  of  this  species  from  the  model  results  in 
a  poorer  fit  to  the  observed  data.  Furthermore,  the 

values  of  and  6m„T  are  reasonably  consistent  with 

f  O  TML .  AH 

previously  published  data  (62,77)  concerning  the  trimethyl¬ 
lead  complexes  of  certain  carboxylic  acids.  Acids  of  pK^ 
similar  to  pKal  of  this  ligand,  such  as  acetylglycine 

(pK  3.40),  have  quite  similar  formation  constants  of  the 
a 

trimethyllead  complex  (the  is  6.6)  as  well  as  similar 

chemical  shifts  (the  SIMPLEX  ^*^9).  These  similarities 
will  be  discussed  in  more  detail  later  in  this  chapter. 

The  large  value  obtained  for  K^g  is  similar  to  that 
calculated  for  the  formation  of  the  2— mercaptoethanol 
complex  of  trimethyllead  (the  K^g  of  which  is  9.05  (-  0.11) 
x  10^) ;  the  deprotonated  sulfhydryl  group  binds  the  tri¬ 


methyllead  cation  in  both  cases.  The  slightly  lower 
magnitude  of  the  N-acetylpenicillamine  binding  constant 
relative  to  that  of  2-mercaptoethanol  is  presumably  due 
to  steric  factors:  the  presence  of  the  bulky  3-methyl 


•  . 
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and  N  acetyl  groups  near  the  sulfhydryl  binding  site  reduces 
the  frequency  of  successful  complexation  with  the  rela¬ 
tively  bulky  trimethyl lead  cation,  compared  to  that  of  the 
small  2 -mercaptoethanol  molecule.  Other  factors,  such  as 
the  larger  pK^  of  the  N-acetylpenicillamine  sulfhydryl 
group  as  well  as  the  additional  electrostatic  attraction 
of  the  deprotonated  carboxyl  group,  tend  to  ameliorate 
the  binding  of  trimethyllead  somewhat.  These  effects 
will  be  discussed  in  more  detail  in  a  later  section. 

4 .  D,L-Penicillamine  Complexes 

(a)  The  Macroscopic  Acid  Dissociation  Constants  of 
D ,L-Penicillamine .  Fully  protonated  D , L-penicillamine 
has  the  structural  formula  shown  by  II: 


hInchcooh 


II 


c(ch3) 2sh 


This  structure  will  be  symbolized  by  H3P+  or  by  an  abbre¬ 


viated  structural  formula  which  portrays  only  the  active 
binding  sites. 

Increasing  the  pH  of  a  solution  containing  H3P  first 
causes  deprotonation  at  the  carboxyl  group.  The  mixed 
acid  dissociation  constant,  K^,  quantifies  the  extent 
of  this  deprotonation  at  a  given  pH. 


r-COOH  K  .  rCOO 

Aal 

-SH  t 


+ 


+  H 


(70) 
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[H2PlaH+ 

[H3p+] 


(71) 


The  deprotonation  of  the  carboxyl  group  is  complete 
well  before  the  deprotonation  of  the  sulfhydryl  and  amino 
groups  occurs.  However,  the  sulfhydryl  and  amino  de¬ 
protonations  are  not  similarly  distinct:  simultaneous 
proton  dissociation  occurs  at  these  sites  over  a  significant 
pH  range.  It  is  not  possible  to  determine  the  individual 
acid  dissociation  constants  of  such  coincident  deprotona¬ 
tions  by  normal  pH  titration  techniques.  Nevertheless, 
the  pH  titration  curve  of  penicillamine  exhibits  separate 
inflection  points  corresponding  to  the  titration  of  two 
acidic  protons  in  this  pH  region.  A  portion  of  both 
acidic  groups  are  titrated  at  the  first  inflection,  resul¬ 
ting  in  stronger  binding  of  the  remaining  proton.  The 
final  proton,  irrespective  of  location,  is  removed  at 
significantly  higher  pH.  The  positions  of  these  inflection 
points  are  used  to  determine  the  macroscopic  acid  dissoci¬ 
ation  constants  K  0  and  K  ^  (97),  but  only  the  overall 

a  2  a  j 

degree  of  protonation  is  considered,  with  no  regard  for  the 
states  of  the  individual  deprotonation  sites. 
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Monoprotonated  (site  unspecified)  and  fully  deprotonated 
penicillamine  are  represented  by  HP~  and  P~~,  respectively. 

In  the  chemical  literature  of  the  last  two  decades, 
the  acid-base  chemistry  of  D ,L-penicillamine  is  usually 
characterized  in  terms  of  the  macroscopic  ionization 
constants  obtained  from  pH  data.  A  survey  of  the  literature 
located  a  number  of  articles  in  which  values  for  these 
constants  had  been  reported.  Several  of  these  were  re¬ 
jected  due  to  inadequate  experimental  descriptions, 
suspected  technical  errors  or  the  employment  of  temperatures 
other  than  25°C  for  the  determinations.  The  remainder  are 
presented  in  Table  II. 

All  of  the  listed  values  were  determined  by  pH 
titration  at  25°C,  but  at  different  ionic  strengths.  This 
invalidates  direct  comparison  of  the  various  ionization 
constants  reported  in  the  literature,  because  the 
values  are  dependent  on  the  ionic  strength.  A  set  of 
corresponding  constants,  corrected  to  an  ionic  strength 
of  0.3  M  (as  employed  in  the  experiments  described  in 
this  thesis),  are  also  presented  in  Table  II.  This  cor¬ 
rection  was  done  using  the  modified  Davies  Equation  (109) 
described  in  King  (110)  and  Robinson  and  Stokes  (111)  : 

1/2 


-log  y .  = 


A  Z2Z2  ^ 


1  +  y 


1/2 


0.1  z1z2u 


(76) 
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The  value  of  A  at  25°C  was  taken  to  be  0.5115  from  Bates 

(112)  . 

Activity  coefficients  were  calculated  using  Equation 
76  for  each  charged  species  in  the  expression  for  the  macro¬ 
scopic  ionization  constants,  both  at  the  reported  ionic 
strength  and  at  0.3  M.  This  determination  of  the  various 
activity  coefficients  enabled  the  conversion  of  the 
reported  ionization  constants  to  the  corresponding  values 
of  those  constants  at  0.3  M  ionic  strength.  An  example 
of  this  procedure  is  outlined  below. 

K  -j ,  as  expressed  in  Equation  75,  is  a  mixed  concen- 
a  j 

tration-activity  ionization  constant.  The  corresponding 

T 

thermodynamic  ionization  constant  K^,  is  defined  by 
Equation  77: 


ap=aH+ 


aHP" 


(77) 


Appropriate  substitution  of  the  general  definition  of 
activity.  Equation  78,  into  the  expression  for  the  thermo¬ 
dynamic  ionization  constant  results  in  Equation  79: 


ai  = 

T  ^p=tp  ]aH+ 
a3  YHp- [HP  ] 


(78) 

(79) 


Recognizing  the  right  hand  side  of  Equation  75  as  a  factor 
of  Equation  79,  Equation  80  can  be  written: 


. 
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Yp= 

yhp" 


(80) 


T 

Since  a  constant,  independent  of  the  ionic  strength 

of  experimental  solutions,  it  follows  that  Equation  81 
expresses  correctly  the  relationship  between  values  of 
Ka3  ot)tained  at  different  ionic  strengths: 


YP~U)  K  =  yP=(2) 

V  *3(1)  Ka3(2) 

'HP  (1)  'HP- (2) 


(81) 


The  bracketted  subscripts  1  and  2  indicate  the  values 
determined  at  the  first  and  second  ionic  strengths,  respec¬ 
tively.  As  stated  earlier,  the  activity  coefficients  at 
both  ionic  strengths  were  calculated  with  Equation  76. 

Substitution  of  the  literature  value  of  K  ,  (at  the  first 

a3 

ionic  strength)  for  allowed  the  calculation  of 

Ka3(2)  value  corrected  to  an  ionic  strength  of  0.3  M)  . 

Similar  calculations  were  performed  to  obtain  corrected 
values  of  K ^  and  Ka2;  the  results  are  presented  in 
Table  II  as  pK  values. 

3. 

Only  two  values  for  pKfll  were  found  in  the  literature, 
together  with  an  unpublished  value  (117)  obtained  using 
the  pmr  chemical  shift  vs.  pH  technique.  The  pK^  value 
of  1.94  was  determined  to  be  the  better  for  several 
reasons:  it  was  the  nearest  to  the  value  of  2.0  obtained 
by  the  pmr  technique,  the  paper  provided  a  much  greater 
quantity  of  experimental  detail  and  the  publication  was 
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the  more  recent  by  a  full  decade.  The  best  values  of 
pK^2  and  P^a3  were  chosen  by  selecting  the  median  of 
each  set  of  corrected  values.  The  median  values  so  ob¬ 
tained  were  the  pK^2  an<^  P^a3  corrected  from  the  results 
published  by  Wilson  and  Martin  (107) . 

(b)  The  Microscopic  Acid  Dissociation  Constants  of 

D , L-Penicillamine .  The  ionizations  of  the  sulfhydryl  and 

amino  groups  of  penicillamine  occur  simultaneously;  the 

macroscopic  constants  K  0  and  K  described  in  the  pre- 

a  z  dj 

vious  subsection  are  composites  of  the  microscopic  ioniza¬ 
tion  constants.  These  microscopic  constants,  denoted  by 
lower  case  "k",  quantify  the  extent  of  deprotonation  of 
the  individual  sulfhydryl  and  amino  groups  as  described 
in  the  accompanying  scheme.  The  subscripts  1,  2  and  3 
refer  to  the  deprotonated  carboxylic  acid,  sulfhydryl  and 
amino  groups,  respectively. 


i-COO 


rCOOH 


K 


al 


-SH 

LnH3 


lnh2 

(P=) 


(82) 
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The  microconstants  defined  by  the  above  scheme  are 
related  to  the  macroconstants  determined  by  pH  titration 


Ka2  k12  +  k13 


K  ^  =  k  ^  +  k 
a3  k123  +  k132 


Rearrangement  of  Equation  84  gives  Equation  85 


Ka3  k123*k132/(k123  +  k132} 


(83) 

(84) 


(85) 


Due  to  the  properties  of  a  cyclic  system,  we  also  have 


Ka2Ka3  k12k123  k13k132 


(86) 


The  derivations  of  these  equations  are  given  in  Edsall 
and  Wyman  (118) . 

Only  three  independent  equations  relate  to  four  micro¬ 
constants;  the  value  of  at  least  one  microconstant  is 
required  to  solve  these  equations  for  the  remaining  micro¬ 
constants.  Using  a  spectrophotometric  method,  Wilson 
and  Martin  (107)  determined  a  value  for  pk-^  °f  8.05  at 
2 5 °C  and  0.16  M  ionic  strength.  This  was  corrected  to 
8.03  at  2 5 °C  and  0.3  M  ionic  strength  with  the  method  of 
calculation  described  in  the  previous  subsection.  Equations 
83,  85  and  86  were  then  employed  to  determine  the  remaining 
microconstants.  At  25°C  and  0.3  M  ionic  strength,  the 
pk  values  of  penicillamine  were  calculated  to  be  8.03 
for  pk^/  8.61  for  pk^,  10.29  for  pk.^3  an<^  9.70  for  Pk^32* 


and  9.70  for  pk 
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These  constants,  together  with  the  pKal  value  of  1.94 
obtained  earlier,  completely  characterize  the  acid-base 
behaviour  displayed  by  D ,L-penicillamine  under  the  experi¬ 
mental  conditions  employed  in  the  study  detailed  by  the 
next  subsection. 

(c)  The  D , L-Penicillamine  Complexes  of  Trimethyllead . 
The  binding  of  trimethyllead  by  D , L-penicillamine  was 
investigated  by  observing  the  chemical  shift  of  the  tri¬ 
methyllead  protons  as  a  function  of  pH.  As  before,  two¬ 
fold  and  four-fold  excesses  of  ligand  were  employed  to 
ensure  accurate  determinations  of  the  constants  of  the 
strong  sulfhydryl  complexes.  The  results  of  these  experi¬ 
ments  are  presented  in  Figure  12. 

In  addition  to  the  trimethyllead  resonances,  signals 
were  observed  for  the  methyl  protons  of  D , L-penicillamine . 
Separate  resonances  were  detected  for  each  methyl  group, 
due  to  the  asymmetry  of  the  neighboring  a-carbon  atom. 

At  any  instant,  the  two  methyl  groups  are  exposed  to 
dissimilar  molecular  environments,  causing  their  resonances 
to  occur  at  different  chemical  shifts  regardless  of  the 
rate  of  rotation  about  the  hindered  bond.  The 

variation  of  the  resultant  chemical  shifts  with  pH  is 
nearly  identical  for  the  1  to  2  and  1  to  4  metal  to  ligand 
solutions;  hence,  only  the  former  is  shown  in  Figure  13. 

The  methyl  resonances,  though  well  separated  at  all  pH 
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pH 


Figure  12.  pH  dependence  of  the  chemical  shift  of  tri 
methyllead  at  25 °C  and  0.3  M  ionic  strength  (maintained 
with  NaC104)  in  solutions  of  0.0100  M  trimethyllead  only 
(A),  0.00498  M  trimethyllead  and  0.00996  M  D/L-penicil- 
1 amine  (  +  )  ,  and  0. 00500  M  trimethyllead  and  0.0199  M 
D,L-penicillamine  (•)• 
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Figure  13.  pH  dependence  of  the  chemical  shift  of  the 
methyl  protons  of  D / L-penicillamine  at  25°C  and  0.3  M  ionic 
strength  (maintained  with  NaClO^)  in  a  solution  of  0.00498 
M  trimethyl lead  and  0.00996  M  D , L-penicillamine . 
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levels,  were  rather  broad  and  low  in  intensity.  There¬ 
fore,  the  measurement  errors  inherent  in  this  data  set 
were  much  higher  than  those  of  the  sharp,  intense  tri- 
methyllead  resonances.  In  addition,  excess  ligand  was 
employed  in  all  these  experiments,  resulting  in  a  smaller 
proportion  of  the  observed  ligand  chemical  shift  being 
due  to  the  shift  of  metal-ligand  complexes.  For  these 
reasons,  the  ligand  resonance  data  was  poorly  suited  for 
the  accurate  calculation  of  formation  constants.  The 
marked  changes  in  chemical  shift  seen  in  Figure  13  are 
for  the  most  part  due  to  deprotonation  of  the  functional 
groups  of  penicillamine;  the  effects  of  complexation  by 
trimethyllead  are  much  less  distinct. 

The  best  fit  to  the  chemical  shift  data  for  solutions 
containing  trimethyllead  and  penicillamine  was  obtained 
with  the  model  which  considered  complexation  at  the 
sulfhydryl  group  only.  An  expanded  model  which  also 
included  a  carboxyl  complex  of  trimethyllead  resulted  in  a 
poorer  fit.  All  models  which  included  2:1  trimethyllead: 
ligand  complexes  resulted  in  no  positive  real  values  for 
the  associated  constants.  This  was  at  least  partially 
attributable  to  the  experimental  conditions,  which  were 
such  as  to  minimize  the  formation  of  (TML)  ^YP6  complexes. 
In  addition,  space-filling  models  indicate  that  the 
formation  of  such  2:1  complexes  with  the  penicillamine 
sulfhydryl  group  is  very  unlikely  due  to  stenc  inter- 


. 
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actions  of  the  trimethyllead  cations  with  the  ligand 
3-methyl  groups.  On  the  other  hand,  trimethyllead-carboxyl 
complexes  are  within  the  realms  of  possibility,  but  are 
not  observed  in  these  data.  This  is  likely  a  result  of 
the  low  basicity  of  the  ligand  carboxylic  acid  group;  a 
correlation  between  pK  and  the  carboxylate  complex  for- 
mation  constant  has  been  observed  for  trimethyllead  (62,77). 
These  points  will  be  discussed  later  in  this  chapter. 

The  most  successful  model  is  depicted  by  Equation  87. 


1 

2 

3 


r-COOH 

-SH 

~NH3 
<H3P  +  ) 


kflS 


(HP‘2) 


11 


13S 


rCOO 


-S-Pb  (CH3) 3 


l-NH 


+ 


132 

pCOO" 

-s” 

Uh2 

(P-) 

y| 

/ 

rCOO 


fl3S 


-S-Pb (CH3) 3 
lNH„ 


(TML.P  ) 


(TML.PH) 


(87) 


The  subscripts  1,  2  and  3  indicate  that  the  product  is 
deprotonated  and  uncomplexed  at  the  carboxyl,  sulfhydryl 
and  amino  groups,  respectively.  A  subscript  f  indicates 
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a  formation  constant  and  a  subscript  S  signifies  that  a 
trimethyl lead  cation  is  bound  to  the  sulfhydryl  group  of 
the  product .  The  formation  constants  k  and  k4ri  are 

I  lb  II  -L  J Id 

thus  defined  as  follows: 


kflS 


[TML.PH] 
[HP~2 ] [TML+] 


kfl3S 


[TML.P  ] 
[P_] [TML+] 


(88) 


(89) 


The  model  for  the  observed  chemical  shift  of  penicillamine 
was  derived  using  Equation  90: 

^OBS  “  ptml6tml  +  ptmloh6tmloh  +  ptmld6tmld 

+  PTML.PH6TML.PH  +  PTML.P-6TML.P"  (90) 

This  model  and  the  appropriate  data  were  supplied  to 
the  curve-fitting  computer  program  KINET  (94),  which 

optimized  the  values  of  kflg,  6TML.PH'  kfi3S  and  6TML.P"‘ 
The  procedure  employed  here  was  quite  analogous  to  that 
described  for  the  metal-ligand  systems  discussed  previously 
in  this  chapter.  These  calculations  yielded  the  following 
optimal  values:  1.13  (-  0.05)  x  10^  for  1.414 

(±  0.001)  for  6tml  pH,  4.22  (±  0.32)  x  105  for  kfl3g  and 
1.373  (±  0.001)  for  6mMT  nwwhere  the  chemical  shifts  are 
in  units  of  ppm  vs.  DSS.  As  before,  the  uncertainties 
are  linear  estimates  of  the  standard  deviation,  calculated 
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by  KINET.  These  results  compare  quite  well  with  the 
other  determinations  in  this  chapter,  given  the  differences 
in  steric  and  electrostatic  interactions  of  the  various 
ligands  with  trimethyllead  cations.  Correlation  of  these 
effects  will  be  discussed  later  in  this  chapter. 

The  value  of  k^g'  the  ionization  constant  of  the 
trimethyllead-complexed  species  TML.PH,  was  calculated 
using  the  relationship: 


:flS  ’  k13S  k123  *  kf  13S 


(91) 


which  is  analogous  to  Equation  86  and  is  derived  in  a 

similar  manner  (118) .  However,  it  must  be  realized  that 

Equation  86  is  defined  for  dissociation  constants  only; 

hence,  the  inverse  of  each  of  the  formation  constants, 

kf IS  and  kfj3S'  is  used  in  the  derivation  of  Equation  91 

-9 

The  value  of  k^3s  was  found  to  be  1.93  x  10  ,  giving 

a  pk  value  of  8.71.  This  is  considerably  less  than 
the  value  of  10.29  determined  for  Pk^23'  but  quite 
similar  to  the  pk13  of  8.61.  Thus  it  may  be  inferred 
that  a  trimethyllead  cation  complexed  to  the  sulfhydryl 
group  of  a  penicillamine  molecule  closely  mimics  a 
similarly  bonded  hydrogen  ion,  in  terms  of  their 
effects  on  the  acid  strengths  of  neighboring  functional 
groups.  The  electrostatic,  inductive  and  steric  factors 
involved  will  be  discussed  later  in  this  work. 
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5 .  L-Cysteine  Complexes 

(a)  The  Macroscopic  Acid  Dissociation  Constants  of 
L~Cy steine .  At  low  pH,  L-cysteine  assumes  the  fully 
protonated  form  shown  by  Ills 

H^NCHCOOH 

I  HI 

CH2SH 

hereafter  symbolized  by  H^C"*"  or  by  an  abbreviated  struc¬ 
tural  formula  identical  to  that  used  for  D,L-penicillamine, 
since  the  functional  groups  are  similarly  situated.  The 
structure  of  III  differs  only  in  the  substitution  of 
hydrogen  atoms  for  the  methyl  groups  adjacent  to  the  sulf- 
hydryl  in  structure  II.  This  apparently  minor  dissimila¬ 
rity  leads  to  significant  differences  in  the  acid-base 
and  complexation  behaviour  of  these  compounds. 

The  macroscopic  ionization  constants  of  L-cysteine 
are  defined  by  the  following  reaction  scheme  and  equations: 


[COOH 
SH 

NH3 

(h3c+) 


rC0°  Ka2 

-SH  «= - 

*-NH3 

(h2o 


r-COO 

-SH 

Uh2 


(HC  ) 


(C  ) 
(92) 


K 


al 


[H2C]aH+ 

[h3c+] 

[HC"]aH+ 

[h2c] 


(93) 


(94) 
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[C  ]  a 


H+ 


[HC  ] 


(95) 


where  [HC  ]  is  the  sum  of  the  concentrations  of  the  two 

possible  monoprotonated  species  shown  in  Equation  92. 

The  rationale  for  this  treatment  of  the  macroscopic 

acid-base  chemistry  of  cysteine  is  identical  to  that  of 

the  previous  section  regarding  penicillamine.  Using  the 

procedure  detailed  therein,  a  series  of  values  for  K 

Ka2'  and  Ka3  were  collected  from  the  literature  and  the 

corresponding  values  of  these  constants  at  0.3  M  ionic 

strength  were  calculated.  The  original  constants  and  those 

corrected  to  an  ionic  strength  of  0.3  M  are  presented  in 

Table  III.  The  median  of  each  set  of  corrected  ionization 

constants  was  used  in  the  calculations  described  later 

in  this  section.  The  most  recent  and  trustworthy  value  of 

_2 

Kal  (113)  was  found  to  be  1.17  x  10  ,  when  corrected  to 

an  ionic  strength  of  0.3  M.  The  median  values  of  an<3 

K  0  were  6.03  x  10”9  and  5.62  x  10”11,  both  of  these 
a3 

calculated  from  the  published  values  of  Coates,  Marsden  and 
Rigg  (119). 


(b)  The  Microscopic  Acid  Dissociation  Constants  of 
L-Cysteine.  The  microscopic  ionization  scheme  of  L-cysteine 
has  been  determined  (119,122)  to  be  similar  to  that  of 
D, L-penicillamine ,  discussed  in  the  previous  section. 
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Equations  83  through  86  thus  apply  for  the  L-cysteine  system 
of  equilibria  as  well  as  for  that  of  D ,L-penicillamine . 

As  before,  these  equations  cannot  be  solved  without  addi¬ 
tional  information. 

Coates,  Marsden  and  Rigg  (122)  obtained  the  required 
data  using  a  spectrophotometric  method,  and  calculated  a 
value  of  8.53  ±  0.01  for  pk^2  at  25°C.  This  corresponds 
to  a  pk-^2  8.38  at  0.3  M  ionic  strength  and  25°C. 

Based  on  this  result,  values  of  8.73  for  pk^,  10.09  for 
pkf2 3  and  9.74  for  pk132  were  computed.  The  previously 
determined  value  of  1.93  for  pK^  completes  the  quanti¬ 
fication  of  the  acid-base  equilibria  described  by  Equation 
96. 

A  comparison  of  the  cysteine  and  penicillamine  micro¬ 
constants  indicates  that  the  sulfhydryl  group  of  the  latter 
is  significantly  more  acidic  than  that  of  cysteine.  The 
ratio  of  the  monoprotonated  penicillamine  species, 

[HP~2] / [HP~31 ,  is  thus  greater  than  the  corresponding 
cysteine  ratio,  [HC-^  1 / [HC^3 ]  .  The  calculated  values  of 


these  ratios  are  3.8  and  2.2,  respectively. 


(c)  The  L-Cysteine  Complexes  of  Trimethyllead . 

As  in  the  earlier  studies,  the  complexation  of  trimethyl¬ 
lead  by  L-cysteine  was  studied  by  observing  the  pH  depen¬ 
dence  of  the  proton  chemical  shift  of  trimethyllead  in 
appropriate  solutions.  Two-fold  and  four-fold  excesses 
of  ligand  were  employed  to  increase  the  sensitivity  of 
the  method  to  the  formation  of  strong  sulfhydryl  complexes, 
at  the  expense  of  the  less  important,  weaker  complexes 
occurring  at  other  functional  groups.  The  results  of 
these  experiments  are  presented  in  Figure  14. 

The  model  which  best  fits  the  observed  data  was  de¬ 
rived  from  a  scheme  of  equilibria  identical  to  that  of 
penicillamine  (Equation  87) ,  except  for  the  inclusion  of 
a  carboxyl  complex  of  trimethyllead. 
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Figure  14 .  pH  dependence  of  the  chemical  shift  of  tri 
methyllead  at  25°C  and  0.3  M  ionic  strength  (maintained 
with  NaC104)  in  solutions  of  0.0100  M  tr imethyllead  only 
(A),  0.00499  M  trimethyllead  and  0.00949  M  L-cysteine 
(4)  and  0.00502  M  trimethyllead  and  0.0190  M  L--cysteine 


(  •)  . 
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All  symbols  are  defined  as  in  the  penicillamine  section 
of  this  chapter,  except  that  the  "P"  of  penicillamine  is 
replaced  with  the  "C"  of  cysteine.  In  addition,  the  sub¬ 
script  "0"  in  indicates  that  the  trimethyllead  cation 

is  bound  to  the  oxygen  atom  of  the  deprotonated  carboxyl 
group . 

In  contrast  to  the  penicillamine  case,  space-filling 
models  indicated  that  an  additional  species,  having  two 
trimethyllead  cations  bound  to  the  cysteine  sulfhydryl, 
was  not  prevented  by  steric  interactions.  However,  while 
a  doubly  coordinated  sulfhydryl  group  is  sterically  possible, 
2:1  trimethyllead-cysteine  complex  formation  is  not  favored 
at  the  low  metal-to-ligand  ratios  used  in  this  work.  In¬ 
corporation  of  this  2:1  species  into  Equation  97  resulted 
in  untenable  values  for  the  corresponding  formation  constant 
and  chemical  shift,  indicating  that  the  contribution  of 
this  species  to  the  observed  chemical  shift  was  too  low 
to  measure  accurately.  With  higher  metal  to  ligand  ratios, 
such  two  to  one  complexes  could  be  considerably  more 
important.  This  speculation  is  reinforced  by  studies  of 
other  workers  (123)  which  indicate  that  two  or  more  of 
the  analogous  (although  less  bulky)  methylmercury  cations 
can  bond  to  one  sulfhydryl  group.  In  any  case,  under 
the  conditions  chosen  for  this  trimethyllead  study, 

Equation  97  best  represents  the  detectable  species  and 
equilibria  in  solution. 


. 
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The  model  of  the  trimethyllead  chemical  shift  in 
the  system  given  by  Equation  97  is  the  following: 

OBS  TML  TML  +  PTMLOH^TMLOH  +  PTMLD^TMLD 

+  PTML.CH2<5TML.CH5  +  PTrtL.CH<5TML.CH 

+  PTML.C"5TML.C“  *98) 

This  model,  together  with  the  appropriate  data,  was  fitted 

to  the  observed  chemical  shift  curve  using  the  computer 

program  KINET  as  described  previously.  The  optimum 

values  of  the  formation  constants  and  the  associated 

chemical  shifts  are  2.2  ±  2.0  for  K f  ,  1.426  ±  0.089  for 

6TMECH5'  9.67  (±  0.27)  x  104  for  kf ±s ,  1.356  ±  0.0005 

for  6tmdch,  9.24  (±  0.50)  x  105  for  kfl3  and  1.325  ± 

0.0005  for  where  the  chemical  shifts  are  in  ppm 

TMI3L 

vs.  DSS.  The  inaccuracies  are  linear  estimates  of  the 

standard  deviations  as  calculated  by  KINET. 

The  large  error  estimates  for  and  result 

directly  from  the  low  concentrations  of  the  carboxyl 

complex:  such  low  concentrations  contribute  relatively 

little  to  the  averaged  trimethyllead  chemical  shift,  from 

which  and  6m..T  are  calculated.  However,  exclusion 

fO  TML.CH2 

of  TMLCH^  from  the  model  gives  a  much  poorer  fit.  In 
addition,  the  calculated  values  of  snd  6 rp^ j<  CH^ 

correlate  reasonably  well  with  other  carboxylate  complex 
data  discussed  in  this  chapter.  For  these  reasons,  the 
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carboxyl  complex  parameters  were  accepted  in  spite  of  the 
poor  precision  observed.  Better  measurements  may  be 
possible  at  higher  metal  to  ligand  ratios. 

The  value  of  k13g  was  calculated  to  be  7.83  x  10-10 
(P^13S  =  9.11)  using  Equation  91.  This  is  much  closer  to 
the  pk13  value  of  8.73  than  to  the  pk123  of  10.09,  due 
to  the  same  factors  which  govern  the  analogous  penicillamine 
deprotonations.  However,  the  trimethyllead  cation  mimics 
the  behaviour  of  a  proton  much  less  closely  for  cysteine 
than  for  penicillamine,  probably  due  to  the  differences  in 
steric  interactions  and  in  inductive  effects  in  the  two 
ligands . 


6 .  Glutathione  Complexes 

(a)  The  Acid  Dissociation  Constants  of  Glutathione. 
When  fully  protonated,  the  tripeptide  glutathione  (y-L- 
glutamyl-L-cysteinylglycine)  has  the  structural  formula 
given  by  IV: 

O  0  0  0 

II  ll  1  ll 

HOCCHCH0CH0CNHCHCNHCH0COH  IV 

I  +  2  2  I 

NH3  CH2SH 

This  molecule  possesses  two  pairs  of  acidic  groups  which 
ionize  simultaneously:  the  two  carboxylic  acid  groups 
which  ionize  in  the  pH  range  0.5  to  6.0  and  the  sulfhydryl 


and  ammonium  groups  which  ionize  in  the  pH  range  7  to  12. 
The  determination  of  the  microscopic  ionization  constants 
of  all  these  groups  was  accomplished  with  studies  of  the 
nmr  chemical  shifts  of  selected  carbon-bonded  protons  as 
a  function  of  pH  (124) .  The  resultant  microscopic 
ionization  scheme  is  given  by  Figure  15.  The  ionization 
reaction  to  which  a  given  microconstant  refers  is  given 
by  the  subscripts:  the  last  number  in  the  subscript 
denotes  the  group  involved  in  the  ionization  step  under 
consideration,  while  the  preceding  numbers  denote  the 
groups  from  which  protons  have  already  ionized.  The 
L-glutamyl  carboxylic  acid  group  is  denoted  by  the  number 
1,  the  glycyl  carboxylic  acid  group  by  2,  the  sulfhydryl 
group  by  3  and  the  ammonium  group  by  4,  as  indicated 
near  H^G+  in  Figure  15.  However,  in  the  case  of  the 
macroscopic  constants  mentioned  below,  the  subscripts 
refer  only  to  the  total  number  of  protons  removed  from 
the  product  of  the  equilibrium,  regardless  of  the  sites 
of  deprotonation. 

The  macroscopic  ionization  constants  exhibited  by 

glutathione  were  calculated  (124)  using  equations  analogou 

to  Equations  83  through  86  of  the  present  work.  The 

micro-  and  macroconstants  of  glutathione  were  presented 

(124)  as  both  p  K  and  pK  values;  only  the  latter  were 

c  a  a 

used  in  this  thesis  and  are  listed  in  Table  IV. 
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TABLE  IV 


IONIZATION  CONSTANTS  OF  GLUTATHIONE  (124) 


MACROSCOPIC3 ' b ' C 

MICROSCOPIC3 ' 

b 

PKal 

2.15 

pkx 

2.19 

+ 

0.04 

Pk2 

3.22 

+ 

0.04 

PKa2 

3.49 

Pfc12 

3.45 

+ 

0.10 

Pk21 

2.42 

+ 

0.01 

?Ka3 

8.76 

pk123 

8.97 

+ 

0.04 

pk124 

•  9.17 

+ 

0.04 

PKa4 

9.56 

pk1234 

9.35 

+ 

0.10 

pk124  3 

9.08 

+ 

0.02 

(a)  Mixed  constants  calculated  in  terms  of  the  hydrogen 
ion  activities,  using  the  pH  meter  readings  directly, 

(b)  At  25°C  and  0.31  to  0.42  M  ionic  strength. 

(c)  Calculated  from  the  microscopic  ionization  constants. 
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(t>)  Glutathione  Complexes  of  Trimethyllead 

(i)  Proton  NMB  Results 

The  complexation  equilibria  of  trimethyllead  and 
glutathione  were  studied  using  measurements  of  the  chemical 
shift  of  the  trimethyllead  protons  as  a  function  of  pH. 

To  obtain  a  greater  precision  for  the  large  formation 
constants,  two-  and  fourfold  excesses  of  ligand  were 
maintained  in  solution.  The  results  of  three  experiments 
are  presented  in  Figure  16. 

The  model  of  metal-ligand  interaction  based  on  the 
large  number  of  different  protonated  forms  shown  in  Figure 
15  is  exceedingly  complex,  involving  simultaneous  complexa¬ 
tion  at  two  pairs  of  deprotonated  functional  groups. 

It  was  determined  that  this  model  exceeded  the  computational 
capabilities  of  the  curve-fitting  program  KINET  (94) . 

In  addition,  the  above-mentioned  low  metal  to  ligand 
ratios  reduced  the  already  small  contributions  of  the 
various  weak  carboxyl  and  2:1  complexes  to  the  observed 
chemical  shift.  These  considerations  necessitated  the 
simplification  of  the  scheme  presented  in  Figure  15. 

The  replacement  of  the  microionization  scheme  of  the 
carboxylic  acid  groups  with  two  successive  macroioniza¬ 
tions  affects  the  sulfhydryl  results  minimally,  since 
both  carboxylates  were  found  to  complex  trimethyllead 
quite  weakly.  In  addition,  complexation  of  the  amino 
group  by  trimethyllead  was  considered  negligible  for 


« 
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Figure  16.  pH  dependence  of  the  chemical  shift  of  tri- 
methyllead  at  25°C  and  0.3  M  ionic  strength  (maintained  with 
NaClO^)  in  solutions  of  0.0100  M  trimethyllead  only  (A), 
0.00501  M  trimethyllead  and  0.00973  M  glutathione  (#), 
0.00501  M  trimethyllead  and  0.0194  M  glutathione  (♦)  and 
0.0749  M  trimethyllead  and  0.146  M  glutathione  (■). 
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reasons  detailed  earlier  in  this  chapter. 

The  simplified  model  which  resulted  in  the  best  fit 
to  the  observed  data  is  given  by  Figure  17.  All  previous 
model  systems  can  be  considered  as  subsets  of  this  system 
of  equilibria.  In  this  figure,  all  subscripts  are  defined 
as  for  Equation  97  and  Figure  15.  In  addition,  the 
subscript  "D"  indicates  the  formation  of  a  two  to  one 
metal  to  ligand  complex,  with  both  trimethyllead  cations 
binding  to  the  sulfhydryl  group  simultaneously.  The  use 
of  brackets  indicates  that  the  model  makes  no  distinction 
between  the  species  so  enclosed. 

The  model  of  the  chemical  shift  behaviour  of  this 
system  is  given  by  Equation  99: 


5obs  ptml5tml  +  ptmloh6tmloh  +  ptmld5tmld 


+  PTML .  GH^TML .  GHlj  +  PTML .  GH2  6TML .  GH. 


+  PTML.GH-6TML.GH-  +  PTML.G  6TML.G: 


+  ptml2.gh6tml2.gh 


(99) 


where  PmikifT  is  defined  by  Equation  100 

TML2 . GH 


P  ^  =  2 [TML0 . GH] /TMLT 

TML2.GH  2 


(100) 


and  all  other  parameters  are  defined  analogously  to  those 
of  models  described  earlier  in  this  chapter. 

Despite  the  simplifications  outlined  above,  the  model 
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given  in  Figure  17  resulted  in  an  exceptionally  good  fit 
to  the  observed  data.  The  optimized  values  of  the  unknown 
parameters  were  the  following:  4.23 (±  0.26)  x  105  for 
kfl24S'  1,335  -  0*0005  for  <5TML  Q=,  3.47  (±  0.05)  x  105 
for  kf 12S '  1,338  -  0.0005  for  6TML  GR_,  10.8  ±  1.7  for 

Kf io '  1*455  1  0,012  for  6tml.gho'  1:2  1  1,1  for  K—' 


1.433  ±  0.076  for  6 


TML.GH$ 


/  17.4  ±  4.9  for  k 


f  12D 


fO 
and 


1*351  -  0.011  for  /  where  all  the  chemical  shifts 

are  in  ppm  vs.  DSS.  The  inaccuracies  quoted  are  the 
linear  estimates  of  the  standard  deviations  calculated 
by  KINET.  The  large  relative  errors  of  the  formation 
constants  and  chemical  shifts  of  the  weak  complexes 
result  from  the  low  concentration  of  these  complexes 
in  the  experimental  solutions.  As  is  the  case  with  weak 
complexes  of  other  ligands  described  earlier  in  this 
chapter,  the  scarcity  of  such  species  hinders  accurate 
determination  of  the  associated  parameters.  However,  the 
small  relative  errors  computed  for  the  sulfhydryl  complexes 
indicate  that  these  are  not  similarly  afflicted  nor  are 
they  significantly  prejudiced  by  the  large  errors  of 
other  constants. 

The  value  of  k124s,  calculated  using  an  equation 
analogous  to  Equation  91,  was  found  to  be  5.45  x  10 

does  not  differ  significantly 

from  the  value  of  9.35  for  P^234;  this  is  reasonable, 
since  the  spatial  separation  ot  the  ionizing  amino  group 


(pk124S  iS  9-26)‘  pk124S 
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and  the  sulfhydryl  complexation  site  is  relatively  large. 

By  a  similar  method,  K^q  was  calculated  to  be 
approximately  2.9  x  10~3  (pKalQ  equals  2.54).  However, 
this  calculation  involves  and  both  of  which 

are  of  limited  precision,  and  does  not  take  into  account 
the  microionization  scheme  of  the  carboxyl  groups  shown 
in  Figure  15.  Therefore,  the  accuracy  of  this  K^q  value 
is  probably  much  less  than  that  of  k124s  above. 

(ii)  Carbon-13  NMR  Results 

To  corroborate  the  results  obtained  above  with  pmr 
techniques,  the  chemical  shifts  of  the  carbon-13  resonances 
of  glutathione  and  trimethyllead  were  studied  as  functions 
of  pH.  The  experimental  conditions  were  similar  to 
those  employed  for  the  pmr  studies,  but  the  extraction 
of  quantitative  data  was  prevented  by  the  formation  of 
small  quantities  of  a  white  precipitate  above  pH  5.5 
for  solutions  containing  0.0947  M  trimethyllead  and  0.146 
M  glutathione.  This  precipitate  turned  brown  and  then 
black  over  a  period  of  several  days.  Such  phenomena 
have  been  consistently  observed  for  solutions  composed 
of  trimethyllead  and  sulfhydryl- containing  ligands  in 
relatively  high  concentrations  at  high  pH.  In  the 
present  study,  the  degree  of  precipitation  was  small 
enough  that  qualitative  interpretation  was  still  possible. 

The  results  for  several  of  the  carbon-13  resonances 
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of  glutathione  and  of  trimethyllead  are  presented  in 
Figures  18,  19  and  20.  The  assignments  of  the  gluta¬ 
thione  resonances  are  those  favored  by  several  authors 
(125-127).  The  carbon  atoms  are  identified  by  Glu,  Cys 
or  Gly,  to  indicate  the  amino  acid  residue  in  which  the 
carbon  is  located,  and  C  ,  Cfl,  C  ,  CONH  or  COOH,  to 

a  3  Y 

indicate  the  particular  carbon  of  that  residue.  The  gluta¬ 
thione  resonances  omitted  from  Figures  18  and  19  did  not 
differ  significantly  from  the  pH  vs.  cmr  chemical  shift 
curves  reported  by  Fairhurst  (128)  for  pure  glutathione 
in  solution. 

No  similar  data  were  available  for  comparison  with 
Figure  20;  nevertheless,  the  glutathione  cmr  data  pro¬ 
vided  confirming  evidence  for  the  model  depicted  in  Figure 
17.  Above  pH  3,  the  observed  cmr  chemical  shift  curves 
for  the  Cys-C^  and  Cys-C^  atoms  are  well  displaced  from 
those  of  pure  glutathione  solutions,  indicating  strong 
complexation  at  the  sulfhydryl  group  down  to  fairly  low 
pH  values.  This  is  consistent  with  the  large  formation 
constants  of  the  glutathione  sulfhydryl  complexes  quanti¬ 
fied  in  the  previous  section.  The  curves  for  the  Glu-COOH 
and  the  Gly-COOH  carbon  atoms  are  not  displaced  signifi¬ 
cantly  in  the  region  between  pH  2  and  6  where  maximum 
carboxyl  complexation  is  calculated  to  occur  (vide  infra: 
Figure  32) .  (The  deviations  observed  for  the  solutions  of 
highest  pH  coincide  with  precipitate  formation.)  This  is 
consistent  with  the  small  formation  constants  measured 
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pH 


Figure  18.  pH  dependence  of  the  carbon-13  chemical  shift 
of  selected  carbons  of  glutathione  at  25°C  in  a  solution 
(117,128)  of  0.21  M  glutathione  only  (dashed  lines)  and  at 
0.4  M  ionic  strength  (maintained  with  NaClO^)  in  a  solution 
of  0.146  M  glutathione  and  0.0947  M  tr imethyllead  (circles). 
Positive  chemical  shifts  indicate  the  resonances  are 


upfield  of  dioxane. 
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pH 

Figure  19.  pH  dependence  of  the  carbon-13  chemical  shift 
of  selected  carbons  of  glutathione  at  25°C  in  a  solution 
(117,128)  of  0.21  M  glutathione  only  (dashed  lines)  and 
at  0.4  M  ionic  strength  (maintained  with  NaClC>4)  in  a 
solution  of  0.146  M  glutathione  and  0.0947  M  trimethyl— 
lead  (circles).  Positive  chemical  shifts  indicate  the 
resonances  are  upfield  of  dioxane. 


■ 
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Figure  20.  pH  dependence  of  the  carbon-13  chemical  shift 
(vs.  dioxane)  of  0.0947  M  trimethyllead  in  solution  with 
0.146  M  glutathione  at  25°C  and  0.4  M  ionic  strength 
(maintained  with  NaClC>4).  Positive  chemical  shifts 
indicate  the  resonances  are  upfield  of  dioxane. 
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for  the  glutathione  carboxyl  complexes  of  trimethyllead . 

The  assumption  that  the  amine  group  of  glutathione  does 
not  complex  detectably  with  trimethyllead  is  borne  out 
by  the  coincidence  of  the  Glu-C^  curves  of  pure  gluta¬ 
thione  and  trimethyllead-complexed  glutathione.  Similarly, 
the  complexation  of  glutathione  by  trimethyllead  does 
not  measurably  enhance  the  acidity  of  the  amide  nitrogens 
of  glutathione.  Thus,  the  cmr  data  tends  to  confirm  the 
model  postulated  above  for  the  complexation  of  trimethyl¬ 
lead  by  glutathione. 

Also  obtained  in  this  study  was  a  value  for  the 

lead-207 :carbon-13  coupling  constant,  which  was  308.8 

Hz  at  pH  5.00  in  the  solution  described  in  this  subsection. 

207  13 

This  is  in  the  vicinity  of  Pb  -  C  coupling  constants 
reported  for  a  number  of  similar  trialkyllead  compounds 
under  a  variety  of  conditions  (81-83). 


C .  Discussion 

A  variety  of  ligands  were  selected  for  the  present 
study  to  permit  the  observation  of  the  effects  of  various 
ligand  parameters  upon  the  formation  of  trimethyllead 
complexes.  Tables  V  and  VI  provide  summaries  of  the 
formation  constants  and  the  chemical  shifts  of  the  various 
ligand- trimethyl lead  complexes  described  in  this  Chapter, 
along  with  some  corresponding  proton,  inorganic  lead  and 
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(e)  pKa  of  the  ionization:  H2NCH2COOH  t  H2NCH2COO”  +  H  ,  from  reference  118,  page  486. 


128 


methylmercury  data  for  comparison.  A  later  section  will 
present  qualitative  explanations  of  the  trends  and 
apparent  anomalies  in  the  tabulated  data  concerning  tri  — 
methyllead  complexes. 

The  results  presented  in  this  chapter  also  demonstrate 
that  the  carboxyl,  sulfhydryl  and  trimethyllead-complexed 
sulfhydryl  groups  of  the  selected  organic  ligands  are 
all  potential  binding  sites  for  (CH^^b*,  and  that  the 
amount  of  complexation  at  a  particular  site  is  strongly 
pH-dependent . 


1 .  The  pH  Dependence  of  the  Complexation  of 
Trimethyllead  by  Organic  Ligands 


(a)  Conditional  Formation  Constants  as  Functions  of 
pH.  The  pH  dependence  of  complex  formation  may  be  illus¬ 
trated  using  conditional  formation  constants;  these  will 
be  distinguished  from  the  formation  constants  described 
earlier  by  a  subscript  c.  A  conditional  formation  constant, 
K^c  ^efi-ne<3  as  the  pH  dependent  equilibrium  constant 
for  the  reaction: 


[TML]  free  +  'LIGANDlfree  Z  [TML •  LIGAND]  tQtal  (101) 


K 


=  [TML. LIGAND] tQtal 
f'°  [™L1  free' [LIGAND1  free 


(102) 


where  [TML] f ree  and  [LIGAND] free  include  all  free  forms  of 
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trimethyl lead  and  ligand,  respectively,  and  [TML . LIGAND] total 
includes  all  forms  (protonated  and  deprotonated)  of  the 
particular  complex  under  study  (131) . 

Although  a  number  of  the  formation  constants  given 
earlier  in  this  chapter  are  large,  the  corresponding  con¬ 
ditional  formation  constants  can  be  quite  small  due  to 
competing  reactions  with  proton  or  hydroxide  ions.  The 
conditional  formation  constants  of  the  trimethyllead 
complexes  with  the  ligands  discussed  in  this  work  are 
presented  as  functions  of  pH  in  Figures  21  through  26. 

Each  K^c  is  almost  always  somewhat  smaller  than  the  cor¬ 
responding  Kf,  and  its  magnitude  is  strongly  pH-dependent . 

The  changes  in  slope  evident  in  the  tails  of  a  few  of  these 
curves  result  from  additional  protonation  or  deprotonation 
at  other  functional  groups  on  the  ligand.  The  differences 
in  amplitude  of  the  curves  of  corresponding  complexes 
are  functions  of  the  differences  in  the  magnitudes  of  the 
various  complex  formation  constants  as  well  as  of  the 
effects  of  the  competing  reactions,  i.e.,  ligand  protonation  and 
trimethyllead  complexation  with  hydroxide  ions.  These 
differences  will  be  explained  in  a  later  section. 

These  Figures  indicate  that  the  one  to  one  trimethyl¬ 
lead  to  sulf hydryl  complex  is  the  strongest  complex  in 
the  pH  range  of  approximately  6  to  12.  At  high  pH, 
there  is  some  dissociation  of  the  complexes  due  to  hydroxide 
ion  competition  for  the  (CH^J^Ph  cation. 


At  intermediate 


pH 


Figure  21.  Conditional  formation  constants  of  the  2- 
mercaptoethanol  complexes  of  trimethyllead  as  function 
of  pH  at  2 5 °C  and  0.3  M  ionic  strength. 
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Figure  22.  Conditional  formation  constants  of  the  glycine 
complexes  of  trimethyllead  as  functions  of  pH  at  25°C 
and  0.3  M  ionic  strength. 
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Figure  23.  Conditional  formation  constants  of  the 
N-acetyl-D,L-penicillamine  complexes  of  trimethyllead  as 
functions  of  pH  at  25°C  and  0.3  M  ionic  strength. 
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Figure  24.  Conditional  formation  constants  of  the 
D/L-penicillamine  complexes  of  trimethyllead  as  functions 
of  pH  at  2 5 °C  and  0.3  M  ionic  strength. 
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Figure  25.  Conditional  formation  constants  of  the 
L-cysteine  complexes  of  trimethyllead  as  functions  of 
pH  at  25°C  and  0.3  M  ionic  strength. 
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Figure  26.  Conditional  formation  constants  of  the  gluta¬ 
thione  complexes  of  trimethyllead  as  functions  of  pH 
at  25 °C  and  0.3  M  ionic  strength. 
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pH,  the  sulfhydryl  complexes  dissociate  due  to  competi¬ 
tive  protonation  of  the  sulfhydryl  group.  For  2-mercapto- 
ethanol  and  glutathione,  small  amounts  of  the  two  to  one 
trimethyllead  to  sulfhydryl  complexes  were  also  detected 
in  this  region.  As  shown  in  Figures  21  and  26,  the  curves 
of  these  two  to  one  complexes  approach  those  of  the  one 
to  one  complexes.  This  occurs  because  the  ratio  of  free 
(CH^)^Pb+  to  deprotonated  sulfhydryl  is  relatively  high 
at  neutral  pH  levels;  hence,  this  is  the  most  favorable 
region  for  the  formation  of  two  to  one  trimethyllead  to 
ligand  complexes.  Such  complexes  were  not  detected  for 
the  majority  of  the  ligands  studied,  due  to  other  factors 
discussed  later  in  this  chapter. 

At  somewhat  lower  pH,  usually  between  pH  2  and  4, 

the  K «-  values  of  the  various  carboxylate  complexes 

f,c 

exceed  all  others.  This  is  due  to  the  much  higher  acidity 
of  the  carboxylic  acid  groups.  Therefore,  at  low  pH 
the  carboxylate  complexes  of  trimethyllead  will  be  the 
strongest  in  the  solution,  even  though  the  carboxylate 
Kf  values  are  several  orders  of  magnitude  less  than  those 
of  the  sulfhydryl  complexes.  At  very  low  pH,  even  the 
car boxy la te  groups  become  protonated  and  trimethyllead 
complexation  declines  rapidly. 

(b)  Species  Distributions  as  Functions  of  pH.  The 
formation  constants  of  trimethyllead  hydroxide  complexes 
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discussed  in  Chapter  III,  the  pK  values  of  the  various 

cl 

ligand  binding  sites  and  the  complex  formation  constants 
determined  earlier  in  this  Chapter  were  used  to  compute 
the  species  distribution  curves  shown  in  Figures  27  through 
32.  These  curves  illustrate  the  pH  dependence  of  the 
fractional  concentration  of  each  species,  where  each 
fractional  concentration  is  defined  as  the  concentration 
of  the  relevant  metal  or  ligand  moiety  found  in  solution 
as  part  of  a  particular  species,  divided  by  the  total 
analytical  concentration  of  the  moiety.  A  trimethyllead 
to  ligand  ratio  of  one  to  two,  common  to  all  the  binding 
studies  described  in  this  Chapter,  was  employed  in  the 
calculation  of  these  curves.  These  calculations  are  also 
dependent  on  the  actual  concentrations  of  trimethyllead 
and  ligand,  which  were  0.005  M  and  0.010  M,  respectively, 
for  each  case. 

The  upper  plot  in  each  of  these  figures  depicts  the 
fractional  concentrations  vs.  pH  of  all  the  species 
containing  trimethyllead  while  the  lower  plots  refer 
strictly  to  those  species  containing  ligand.  At  low  pH, 
the  primary  species  in  solution  are  (CH^)-^^^  and  fully 
protonated  ligand.  At  higher  pH  levels,  a  succession  of 
complexes  form,  according  to  the  values  of  the  various 
formation  constants  and  pK  's.  The  proportions  found  at 

d 

each  pH  value  are  consistent  with  the  conditional 
formation  constant  data.  In  addition,  these  plots  provide 
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Figure  27.  pH  dependence  of  the  trimethyllead  (upper  plot) 
and  2-mercaptoethanol  (lower  plot)  species  distribution  in 
a  solution  containing  0.005  M  trimethyllead  and  0.010  M 
ligand  at  25°C  and  0.3  M  ionic  strength.  Symbols  are 
defined  as  in  Equations  13  and  40. 
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Figure  28.  pH  dependence  of  the  tr imethyllead  (upper  plot) 
and  glycine  (lower  plot)  species  distribution  in  a  solution 
containing  0.005  M  trimethyllead  and  0.010  M  ligand  at  25°C 
and  0.3  M  ionic  strength.  Symbols  are  defined  as  in 
Equations  13  and  55  -  60. 
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Figure  29.  pH  dependence  of  the  tr imethyllead  (upper  plot) 
and  N-acetyl-D ,L-penicillamine  (lower  plot)  species 
distribution  in  a  solution  containing  0.005  M  trimethyl- 
lead  and  0.010  M  ligand  at  25°C  and  0.3  M  ionic  strength. 
Symbols  are  defined  as  in  Equations  13  and  66. 
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Figure  30.  pH  dependence  of  the  trimethyllead  (upper  plot) 
and  D , L-penicillamine  (lower  plot)  species  distribution 
in  a  solution  containing  0.005  M  trimethyllead  and  0.010  M 
ligand  at  25°C  and  0.3  M  ionic  strength.  Symbols  are 
defined  as  in  Equations  13  and  87. 
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Figure  31.  pH  dependence  of  the  trimethyl lead  (upper  plot) 
and  L-cysteine  (lower  plot)  species  distribution  in  a 
solution  containing  0.005  M  trimethyllead  and  0.010  M 
ligand  at  25°C  and  0.3  M  ionic  strength.  Symbols  are 
defined  as  in  Equations  13  and  97. 
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Figure  32.  pH  dependence  of  the  tr imethyllead  (upper  plot) 
and  glutathione  (lower  plot)  species  distribution  in  a 
solution  containing  0.005  M  trimethyllead  and  0.010  M 
ligand  at  25°C  and  0.3  M  ionic  strength.  Symbols  are 
defined  as  in  Equation  13  and  Figure  17. 


information  regarding  the  degree  of  trimethyllead  complexa 
tion  to  certain  functional  groups  for  differing  ligand 
protonation  states.  Penicillamine,  cysteine  and  gluta¬ 
thione  exhibit  separate  curves  for  the  sulfhydryl  com¬ 
plexes  protonated  and  deprotonated  at  the  amino  group, 
as  does  glycine  for  its  analogous  carboxylate  complexes. 
Glutathione  also  shows  separate  curves  for  the  glycyl 
and  glutamyl  carboxylate  complexes  of  trimethyllead,  with 
the  glycyl  complex  dominating  above  pH  3.  (The  proportion 
of  any  trimethyllead-complexed  species  in  the  lower  plots 
is  limited  to  a  maximum  of  0.5  due  to  the  two-to-one 
excess  of  ligand  over  metal  concentration.) 

The  carboxylate  complex  of  fully  deprotonated  glycine 
with  trimethyllead  is  shown  in  Figure  28  to  be  reduced 
in  proportion  by  competition  for  the  trimethyllead  cation 
with  the  dimeric  species  (  (CH^) ^Pb) 20H  *  T^e  contrast 
in  total  proportions  between  this  complex  and  the  various 
sulfhydryl  complexes  shown  in  the  other  species  distri¬ 
bution  plots  is  quite  striking,  supporting  the  assumption 
made  earlier  in  this  Chapter  that  the  concentration  of 
such  carboxyl  complexes  is  negligible  in  amino  acid 
solutions  containing  an  excess  of  sulfhydryl  groups. 
However,  under  conditions  of  excess  trimethyllead,  detec 
tion  and  quantitation  of  these  fully  deprotonated  carboxyl 
complexes  may  be  possible. 

The  amounts  of  two-to-one  trimethyllead- to -2-mercapto 
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ethanol  or  glutathione  complexes  are  also  seen  to  be 
small,  due  to  competition  with  several  other  equilibria. 
However ,  the  proportion  of  such  complexes  in  solution 
is  highly  dependent  on  the  total  concentration  of  tri- 
methyllead:  at  higher  metal  concentrations  and  at 

appropriate  pH  levels,  these  species  may  become  major 
solution  components.  Under  such  conditions,  similar  two- 
to-one  complexes  may  also  form  in  measurable  quantities 
with  some  of  the  remaining  sulf hydryl-containing  ligands. 

The  composition  of  a  trimethyllead-ligand  solution 
can  be  seen  in  Figures  21  through  32  to  be  strongly  pH 
dependent.  Therefore,  the  control  of  pH  enables  selection 
of  the  type  of  trimethyllead  complex  occurring  in  solution. 
For  example,  at  physiological  pH  values  (near  7.4),  the 
major  species  present  under  conditions  of  excess  ligand 
will  be  the  one-to-one  trimethyllead-sulf hydryl  complex. 

The  two-to-one  and  the  carboxylate  complexes  of  trimethyl¬ 
lead  are  not  likely  to  occur  in  significant  proportions 
under  these  conditions,  unless  most  of  the  free  sulf hydryl 
groups  are  subsequently  rendered  relatively  inactive  by 
some  other  means,  such  as  oxidation. 

2 .  Other  Factors  Influencing  the  Complexation  of 

Trimethyllead  by  Organic  Ligands 

The  trimethyllead  cation  is  a  Lewis  acid;  hence,  the 
strength  of  its  complexes  will  be  dependent  upon  the 
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inherent  basicities  of  the  ligand  binding  sites  involved. 

The  formation  constants  for  all  one-to-one  complexes 
described  in  this  Chapter  are  plotted  versus  this  basicity 
(expressed  in  terms  of  pK  values)  in  Figure  33.  Included 

a 

in  this  plot  are  the  formation  constants  for  several 
trimethyllead-carboxylic  acid  complexes  (62) .  In  this 
Figure,  reasonably  strong  correlation  of  log  with  pK 
is  demonstrated  by  the  various  carboxylate  complexes  and, 
to  a  lesser  extent,  by  the  sulfhydryl  complexes  as  well. 
There  are,  nonetheless,  a  number  of  factors  including 
differences  in  polarizability,  in  steric  and  conformational 
requirements  and  in  electrostatic  interactions,  which 
cause  significant  deviations  from  linearity  in  the  data. 

Firstly,  the  central  lead (IV)  atoms  of  trimethyllead 
are  much  larger  and  more  polarizable  than  are  protons 
(although  the  polarizability  of  this  lead  species  is  much 
less  than  that  of  inorganic  lead).  As  noted  in  Chapter  I, 
the  bonds  of  trimethyllead  with  polarizable  ligand  atoms, 
such  as  the  sulfhydryl  sulfur,  will  have  significant 
covalent  character  and  may  include  pH-dll  bonding  as  well. 

On  the  other  hand,  the  weak  bonds  of  trimethyllead  with 
ligand  atoms  of  low  polarizability  ("hard  Lewis  bases, 
in  the  sense  described  by  Pearson  (64))  will  have  little 
covalent  character  and  no  p IT — d IT  bonding.  Hence,  the  bonding 
in  the  weak  carboxylate  complexes  of  trimethyllead  is 
primarily  ionic.  The  proton  is  a  much  harder  Lewis  acid 
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Figure  33.  The  dependence  of  the  logarithm  of  the  tri- 
raethyllead  complex  formation  constant  upon  the  pK&  values 
of  the  respective  free  ligands.  Complexes  shown  are  those 
studied  in  this  work  (O)  and  in  references  62  (0)  and 


63  (□)  . 
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than  trimethyllead ,  thus  forming  stronger  ionic  bonds 
with  hard  bases  than  does  trimethyllead  (64).  Interesting¬ 
ly,  the  "soft"  methylmercury,  the  "borderline"  Pb++  and 
the  hard  hydrogen  ions  all  form  stronger  bonds  with 
sulfhydryl  and  carboxyl  groups  than  does  trimethyllead, 
indicating  the  importance  of  factors  other  than  hard  and 
soft  Lewis  acidity  and  basicity  in  the  formation  of  tri¬ 
methyllead  complexes.  (However,  it  should  be  noted  that 

2+ 

the  strengths  of  the  Pb  and  trimethyllead  complexes 
cannot  be  compared  directly  because  the  ligands  studied 
in  this  work  may  chelate  Pb  using  several  donor  groups 
simultaneously,  whereas  trimethyllead  reacts  as  a  one- 
coordinate  species  only.) 

Steric  factors  have  a  sizeable  influence  on  the 
complexation  behaviour  of  trimethyllead.  The  presence 
of  the  three  equatorial  methyl  groups  severely  reduces 
the  probability  of  the  formation  of  a  chemical  bond  upon 
collision  of  trimethyllead  and  ligand  molecules.  Protons, 
lead (II)  ions  and  even  methylmercury  cations  are  much 
less  hindered  than  trimethyllead  cations.  Hence,  it  is 
not  surprising  that  these  form  stronger  complexes  with 
the  ligands  of  this  study  than  does  trimethyllead.  This 
effect  is  amplified  by  the  presence  of  bulky  non-reactive 
groups  on  some  ligand  molecules.  For  example,  molecular 
models  indicate  that  the  3-niethyl  groups  of  penicillamine 
and  N-acetylpenicillamine  interact  sterically  with  the 
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trimethyl lead  methyl  groups  such  that  the  complex  formation 
constants  are  significantly  lower  than  those  of  the  analo¬ 
gous  cysteine  complexes.  The  N-acetyl  group  also  hinders 
complex  formation  to  some  extent. 

A  complicating  factor  is  that  of  electrostatic 
attraction  and  repulsion  between  the  .trimethyllead  cation 
apd  the  cationic  or  anionic  ligand  groups  remote  from 
the  site  of  complexation.  Such  interactions  are  generally 
small,  but  are  nonetheless  noticeable  if  the  interfering 
ionic  group  is  situated  quite  near  to  the  complexing 
functional  group.  This  is  especially  true  if  steric 
factors  increase  the  probability  of  these  interactions. 

For  example,  the  formation  of  the  monoprotonated  (protonated 
at  the  ammonium  group)  penicillamine  complex  of  trimethyl¬ 
lead  is  significantly  reduced  due  to  a  combination  of 
electrostatic  and  steric  effects.  Space-filling  model 
studies  indicate  that  the  conformations  allowing  closest 
approach  of  trimethyllead  to  the  ligand  are  those  in 
which  the  ammonium  group  is  situated  near  the  sulfhydryl 
binding  site,  thus  optimizing  the  repulsive  electrostatic 
interaction.  Conversely,  the  monoprotonated  cysteine 
molecule  may  adopt  conformations  minimizing  electrostatic 
repulsion  of  cationic  trimethyllead;  the  formation  constant 
of  the  monoprotonated  cysteine  complex  is  much  larger. 

Another  example  of  electrostatic  effects  is  that 
of  the  monoprotonated  glutathione  complex.  In  this  case, 
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the  glutamyl  ammonium  group  is  well  removed  from  the 
sulfhydryl  binding  site,  resulting  in  no  significant  reduc¬ 
tion  of  the  sulf hydryl-trimethyllead  complex  formation 
constant.  Similarly,  the  glutathione  glycyl  carboxylate 
complex  of  trimethyllead  is  well  removed  and  immune  to 
the  effects  of  the  ammonium  group.  However,  the  glutamyl 
carboxylate  complex  forms  adjacent  to  the  glutamyl 
ammonium  group:  this  complex  forms  only  to  a  very  small 
extent . 

A  further  example  of  the  effect  of  a  nearby  positively 
charged  protonated  group  is  provided  by  the  glycine  system. 
The  two  complexes  most  probably  involved  are  (CH^^Pb- 
C^CC^NH^  and  (CH^ )  2^t)“02CCH2NH2 .  The  difference  observed 
between  the  values  of  the  two  formation  constants  is  due 
to  the  protonation  state  of  the  amino  group. 

Other  examples  of  electrostatic  interactions  include 
the  penicillamine  carboxylate  and  the  two  N-acetyl- 
penicillamine  complexes  of  trimethyllead.  The  concen¬ 
tration  of  the  penicillamine  complex  is  below  experimental 
detection  limits  due  to  steric  hindrances  and  the  afore¬ 
mentioned  sterically-optimized  electrostatic  repulsion  of 
cationic  trimethyllead.  The  N-acetyl-penicillamine 
complexes  are  free  of  electrostatic  repulsion,  but  the 
acetyl  blocking  group  causes  increased  steric  interactions 
between  the  ligand  and  trimethyllead. 

Steric  and  electrostatic  interactions  also  affect 
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the  formation  of  two— to— one  trimethyllead— to —ligand 
complexes.  Only  those  complexes  which  are  relatively  free 
of  such  interferences  at  the  sulfhydryl  binding  site  will 
form  measurable  quantities  of  such  complexes.  In  this 
study ,  the  only  ligands  found  to  be  of  this  type  were 
2— mercaptoethanol  and  glutathione.  No  fully  deprotonated 
two-to-one  complexes  were  observed,  due  to  the  competition 
of  hydroxide  ion  for  the  trimethyllead  cation  at  the 
high  pH  values  required  to  fully  deprotonate  the  ligands 
in  this  study. 

Table  V  presents  overall  formation  constants  for 
some  complexes  of  inorganic  lead(II).  However,  comparison 
of  the  trimethyllead  results  is  somewhat  misleading  without 
a  consideration  of  the  chelation  effect.  Penicillamine 
and  cysteine  probably  act  as  tridentate  ligands  toward 
Pb++  (114,121),  resulting  in  larger  overall  formation 
constants  and  more  complicated  behaviour  than  that  exhi¬ 
bited  by  trimethyllead.  (Glutathione,  however,  appears 
to  be  monodentate  towards  Pb+  (125).)  In  contrast,  it 
has  been  found  in  this  work  and  in  other  studies  (62,63, 
77,129)  that  trimethyllead  reacts  as  a  one-coordinate 
species  with  a  variety  of  ligands;  chelation  is  prevented 
by  the  equatorial  methyl  groups.  Thus,  the  strength  of 
trimethyllead  complexes  is  significantly  less  than  the 
analogous  Pb"*~+  complexes.  In  addition,  the  oxidation 
states  of  the  lead  atoms,  the  overall  charges  and  the 
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stGric  effects  are  different  for  trimethyllead  and  Pb++. 
The  low  degree  of  correlation  observed  between  the  for¬ 
mation  constants  of  analogous  complexes  of  trimethyllead 
and  Pb++  is  therefore  to  be  expected. 

In  contrast,  the  methylmercury  complex  formation 
constants  shown  in  Table  V  are  somewhat  more  easily  cor¬ 
related  to  those  of  trimethyllead.  This  is  because  the 
central  Pb(IV)  of  trimethyllead  is  isoelectronic  with 
the  Hg(II)  of  methylmercury,  the  two  metal  ion  species 
have  identical  +1  charges  and  both  species  bind  to  various 
ligand  functional  groups  in  primarily  one-to-one  ratios. 
However,  steric  factors  are  relatively  unimportant  in 
methylmercury  complexation ,  and  in  addition,  methylmercury 
is  a  very  "soft",  polarizable  Lewis  acid  (64)  compared 
to  trimethyllead.  Therefore,  methylmercury  binds  much 
more  strongly  with  a  wide  variety  of  ligands  (especially 
"soft"  Lewis  bases  such  as  sulfhydryl  groups)  than  does 
trimethyllead,  despite  their  above-mentioned  similarities. 

The  correlation  of  the  Lewis  basicity  of  the  ligand 
functional  group  to  the  nmr  chemical  shift  of  the  tri¬ 
methyllead  protons  in  the  complex  may  be  examined  in  the 
plot  of  6  vs.  pK  shown  by  Figure  34.  The  Lewis 

pj-uu  uj.  0C0MPLEX  f  a 

basicity  of  an  electron  donor  group  is  by  definition  a 
measure  of  its  electron  donating  power:  certain  electrons 
of  the  base  move  towards  an  acidic  center  by  an  amount 
proportional  to  the  strength  of  the  bond.  Increased 
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Figure  34.  The  dependence  of  the  chemical  shift  of  the 
trimethyllead  methyl  protons  in  certain  trimethyllead 
complexes  upon  the  pK  values  of  the  respective  free  ligands. 
Complexes  shown  are  those  studied  in  this  work  (O)  and 
in  references  62  (0)  and  63  (□). 
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bond  strength  results  in  increased  electron  density  near 
the  acidic  center,  which  in  the  case  of  trimethyllead 
results  in  increased  electron  density  near  the  electron- 
withdrawing  methyl  groups.  The  protons  of  these  methyl 
groups  are  therefore  better  shielded  from  the  external 
magnetic  field  in  an  nmr  experiment,  causing  smaller 
chemical  shifts  relative  to  DSS  for  the  stronger  complexes. 
The  ligands  which  form  stronger  complexes  with  trimethyl¬ 
lead  in  general  form  stronger  complexes  with  hydrogen 
ions  as  well,  due  to  the  greater  basicity  of  the  functional 
group. 

A  number  of  deviations  from  linearity  are  observed 
in  Figure  34.  Other  factors,  such  as  steric  interactions, 
which  are  unrelated  to  the  ligand  basicity,  can  nonetheless 
reduce  the  strength  of  the  trimethyllead  complex.  This 
reduces  the  shielding  of  the  trimethyllead  protons  and 
increases  the  chemical  shift  because  the  average  metal- 
ligand  bond  length  is  increased.  These  weaker  complexes 
undergo  a  lesser  degree  of  electron  transfer  from  base  to 
acid.  This  increase  in  chemical  shift  is  observed  for 
the  sterically-hindered  penicillamine  and  N-acetylpenicil- 
1 amine  sulfhydryl  complexes  in  Figure  34.  Electrostatic 

effects  appear  to  be  negligible. 

The  points  presented  in  the  above  discussions  are 
necessarily  qualitative,  but  they  nonetheless  provide 
consistent  explanations  of  the  nature  of  trimethyllead 
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binding  with  the  various  ligands  studied  in  this  thesis. 

3 .  Conclusions 

It  should  be  possible  to  use  the  material  presented 
in  this  thesis  to  elucidate  the  behaviour  of  trimethyllead 
in  biological  systems.  For  example,  at  the  physiological 
pH  level  of  7.4,  the  results  of  this  thesis  when  combined 
with  previous  studies  (62,63,77,129)  indicate  the  preferred 
site  of  trimethyllead  binding  to  be  the  sulfhydryl  group 
of  sulfhydryl-containing  ligands.  The  blocking  of  sulf¬ 
hydryl  active  sites  on  proteins  of  various  sorts  is 
probably  a  major  avenue  of  trimethyllead  poisoning,  as  it 
is  with  methylmercury  (84).  However,  the  evidence  is  not 
conclusive  because  the  body  of  research  on  trimethyllead 
complexation  by  biological  ligands  is  quite  small.  Studies 
on  the  binding  of  trimethyllead  by  imidazole  and  phenolic 
hydroxide  groups,  for  example,  could  also  shed  light  on 
the  mechanisms  of  trimethyllead  toxicity.  A  detailed 
knowledge  of  the  complexation  behaviour  of  trimethyllead 
in  aqueous  and  non-aqueous  systems  will  be  required  to 
develop  a  complete  understanding  of  the  physiological 
effects  of  trimethyllead  intoxication  and  to  optimize  the 
design  of  appropriate  therapeutic  agents. 
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